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Honorable Commissioner of Patents: 
I, Xu Youhao, hereby declare that: 

1. I'm living in Beijing, People's Republic of China. 

2. I'm a citizen of People's Republic of China. 

3. From September 1981 to July 1985, I studied in the Department 
of Petroleum Refining in East China University of Science and 
Technology and obtained my Bachelor's degree; from September 1985 to 
July 1987, I studied in the Department of Chemical Engineering in East 
China University of Science and Technology and obtained my Master's 
degree; from September 2001 to December 2006, I studied in the 
Research Institute of Petroleum Processing and obtained my Doctor's 
degree. 

4. Since August 1987, I've been devoting myself to the development 
of catalytic cracking technology in the Research Institute of Petroleum 
Processing, and have twenty years' experience. I took part in the 
development of Deep Catalytic Cracking (DCC) technology, and I'm 
currently in charge of the development and application of a newly series 
of FCC processes for producing maximum iso^araffins in cracked 
naphtha or clean gasoline and propylene (MIP and MIP-CGP). I have 
written more than forty papers on this subject, six of which were 
published in internationally significant periodicals and international 
academic conferences, and applied for, altogether, more than seventy 
patents domestically and abroad. 

5. In US Patent Application 09/553,990, I propose a novel riser 
reactor, wherein the ratio of the diameter of a second reaction zone to the 
diameter of a first reaction zone ranges from 1.5: 1 to 5: 1. However, in 
catalytic cracking of the prior art, the ratio of the diameter of a second 
reaction zone to the diameter of a first reaction zone in a riser reactor 
should be less than 1.5: 1. 



6. On the basis of the European patent EP0171460, Kmecak et al., 
the Examiner in charge of this apphcation drew the conclusion, after 
manually making a measurement in Figure 8 of EP0171460, that "DR2: 
DRi is equal to approximately 3: 1". In my point of view, this is 
impossible because of the following reasons: 

6.1. In catalytic cracking of the prior art, the ratio of the diameter of 
a so-called second reaction zone to the diameter of a so-called first 
reaction zone in a riser reactor should be less than 1 .5: 1 . 

6.2. Figure 8 in EP0171460 is merely a schematic drawing without 
the specific dimensions available. Thus, the value obtained by manual 
measurement on the basis of the schematic drawing is inaccurate. One 
skilled in the art would never do this. 

6.3. If "DR2: DRi is equal to approximately 3: 1", it will be obtained 
upon scientific calculation (see Item 8 for details) that the ratio of the 
velocity at the outlet of the first reaction zone to the velocity at the inlet 
of the first reaction zone, i.e. U]/uo, is 18.7, which is far greater than the 
generally recognized ratio of less than 3 in the art. Thus, the Examiner's 
conclusion cannot be fulfilled. 

7. In my point of view, the present invention achieves unexpected 
technical effects, which can be seen clearly from Example 1 and the 
corresponding Comparative Example 1, and Example 2 and the 
corresponding Comparative Example 2, as given hereinbelow. 

Comparative Example 2 is a complementary test, and Example 1, 
Comparative Example 1 and Example 2 are recited in the original 
application document. The properties of feedstock used in the Examples 
and Comparative Examples and the properties of the catalyst are 
respectively listed in Tables 1 and 2. 

Example 1 

This Example showed the use of catalyst A in a novel pilot plant 
riser reactor according to the present invention for the production of 
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isobutane and isoparaffin enriched gasoline. 

The total height of the prelift zone, the first reaction zone, the 
second reaction zone, the outlet zone of the riser reactor was 1 5 meters, in 
which the height of the prelift zone with the diameter of 0.025 meter was 
1.5 meters, the height of the first reaction zone with a diameter of 0.025 
meter was 4 meters, the height of the second reaction zone with a 
diameter of 0.1 meter was 6.5 meters, the height of the outlet zone with a 
diameter of 0.025 meter was 3 meters. The isotrapezia vertex angle of the 
vertical section of the conjunct section between the first reaction zone and 
the second reaction zone was 45°. The isotrapezia base angle of the 
vertical section of the conjunct section between the second reaction zone 
and the outlet zone was 60°. 

The reactor and operating conditions were different from those used 
in conventional catalytic cracking, while the regenerating and 
fractionating systems were the same as those used in conventional 
catalytic cracking. The operating conditions and product slate are listed in 
Table 3, and the gasoline properties are listed in Table 4. 

Comparative Example 1 

Compared with example 1, this comparative example was practiced 
in a conventional pilot plant iso-diameter riser reactor. 

The operating conditions and product slate are listed in Table 3, and 
the gasoline properties are listed in Table 4. 

It can be seen clearly from Tables 3 and 4 that: 

1) in Example 1, 35.07 wt% of LPG was isobutane, while in 
Comparative Example 1, 15.74 wt% of LPG was isobutane (see Table 3); 

2) in Example 1, the gasoline had an isoparaffin content of 36.0 wt% 
and an olefin content of 28.11 wt%, while in Comparative Example 1, the 
gasoline had an isoparaffin content of 11.83wt % and an olefin content of 
56.49 wt% (see Table 4). 
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Example 2 



This Example showed the use of gasoline with high olefin content as 
a quenching medium in a novel riser reactor for the production of 
isoparaffin enriched gasoline by means of catalytic conversion in the riser 
reactor. 

The total height of the prelift zone, the first reaction zone, the 
second reaction zone, the outlet zone of the riser reactor was 15 meters, in 
which the height of the prelift zone with the diameter of 0.025 meter was 
1.5 meters, the height of the first reaction zone with a diameter of 0.025 
meter was 4 meters, the height of the second reaction zone with a 
diameter of 0.05 meter was 6.5 meters, the height of the outlet zone with 
a diameter of 0.025 meter was 3 meters. The isotrapezia vertex angle of 
the vertical section of the conjunct section between the first reaction zone 
and the second reaction zone was 45°. The isotrapezia base angle of the 
vertical section of the conjunct section between the second reaction zone 
and the outlet zone was 60°. 

The catalyst and gasoline feedstock used in this Example were the 
same as those used in Example 1 . The gasoline produced in comparative 
example 1 as a quenching medium was charged into the conjunct section 
between the first reaction zone and the second reaction zone. The other 
steps were basically the same as those in Example 1 . 

The operating conditions and product slate are listed in table 5, and 
the gasoline properties are listed in table 6. 

Comparative Example 2 

Compared with example 2, this comparative example was practiced 
in a conventional pilot plant iso-diameter riser reactor. 

The operating conditions and product slate are listed in Table 5, and 
the gasoline properties are listed in Table 6. 

It can be seen clearly from Tables 5 and 6 that, when the ratio of the 
diameter of the second reaction zone to the diameter of the first reaction 
zone was 1 in Comparative Example 2, LPG had an isobutane content of 



only 14.63 wt% and the gasoline had an isoparaffin content of only 9.51 
wt%; when the ratio of the diameter of the second reaction zone to the 
diameter of the first reaction zone was 2 in Example 2, LPG had an 
isobutane content of up to 34.15 wt% and the gasoline had an isoparaffin 
content of up to 43.86 wt%, which are respectively 19.52 percentages and 
34.35 percentages higher than those in Comparative Example 2. 

In other words, when the ratio of the diameter of the second reaction 
zone to the diameter of the first reaction zone is 2 or 4, the present 
invention has a significant advantage in respect of effects over the prior 
art in which the diameter ratio is 1 : 1 . 



Table 1 



Density (20 °C), kgW 


890.5 


Kinematic Viscosity (100 °C), mm^/s 


5.08 


Carbon Residue, wt% 


0.7 


Pour Point, °C 


40 


Total Nitrogen, wt% 


0.16 


Sulfiir, wt% 


0.53 


Carbon, wt% 


85.00 


Hydrogen, wt% 


12.62 


Heavy Metal Content, ppm 




Ni 


0.16 


V 


0.15 


Fe 




Cu 




Na 


0.45 


Distillation, °C 




IBP 


278 


10% 


385 


30% 


442 


50% 


499 


70% 




90% 




EP 
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Table 2 



Catalyst Name 


A 


B 


Trade Mark 


ZCM-7 


CRP-1 


Chemical Composition, wt% 






Aluminum Oxide 


46.4 


54.2 


Sodium Oxide 


0.22 


0.03 


Ferric Oxide 


0.32 




Apparent Bulk Density, kg/m^ 


690 


860 


Pore Volume, ml/g 


0.38 


0.26 


Surface Area, mVg 


164 


160 


Attrition Index, wt% hr"' 




1.2 


Particle Size Distribution, 
wt% 






0~40 microns 


4.8 


26.0 


40—80 microns 


47.9 


60.8 


>80 microns 


47.3 


13.2 



Table 3 





Example 1 


Comparative 
Example 1 


Reactor 


Novel Riser of 


Conventional 




the Present 


Riser 




Invention 




Reaction Temperature, °C 




495 


The First Reaction Zone 


545 


— 


The Second Reaction Zone 


495 


— 


Reaction Time, sec 


5.0 


2.89 


The First Reaction Zone 


1.0 


— 


The Second Reaction Zone 


3.5 


— 


The Outlet Zone 


0.5 


— 


C/O Ratio 


4.5 


4.5 


S/O Ratio 


0.05 


0.05 


Product Slate, wt% 






Dry Gas 


1.83 


1.62 


LPG 


16.11 


11.88 


In which Isobutane 


5.65 


1.87 


Gasoline 


46.86 


41.59 


LCD 


23.44 


22.81 


HCO 


7.77 


18.76 


Coke 


3.88 


2.86 


Loss 


0.11 


0.48 



Table 4 





Example 1 


Comparative 
Example 1 


Reactor 


Novel Riser of 


Conventional 




the Present 


Riser 




Invention 




Density (20 "C), kgW 


743.6 


749.8 


Octane Number 






RON 


90.0 


91.0 


MON 


79.0 


79.8 


Induction Period, min 


>1000 


>485 


Existent Gum, mg/lOOml 


2.0 


2.0 


Sulfur, wt% 


0.0095 


0.0120 


Nitrogen, wt% 


0.0028 


0.0033 


Carbon, wt% 


86.14 


86.81 


Hydrogen, wt% 


13.72 


13.12 


Distillation, °C 






IBP 


46 


50 


10% 


73 


77 


30% 


95 


99 


50% 


114 


122 


70% 


143 


145 


90% 


171 


175 


EP 


202 


205 


Gasoline composition, wt% 






Paraffins 


41.01 


15.81 


In which Iso-paraffins 


36.00 


11.83 


Naphthenes 


7.20 


6.50 


Olefins 


28.11 


56.49 


Aromatics 


23.68 


21.20 



Table 5 





Example 2 


Comparative 
Example 2 


Reactor 


Novel Riser of 


Conventional 




the Present 


Riser 




Invention 




Operating Conditions 






Reaction Temperature, °C 




495 


The first reaction zone 


545 




The second reaction zone 


495 




Reaction Time, sec 


5.3 


2.8 


The first reaction zone 


0.8 




The second reaction zone 


3.9 




The outlet zone 


0.6 




C/O ratio 


5.0 


5.0 


S/O ratio 


0.05 


0.05 


Product Slate, wt% 






Dry Gas 


1.78 


1.65 


LPG 


17.51 


11.07 


In which Iso-butane 


5.98 


1.62 


Gasoline 


47.98 


41.39 


LCO 


22.30 


22.81 


HCO 


6.22 


19.76 


Coke 


4.00 


3.00 


Loss 


0.21 


0.32 



Table 6 





Example 2 


Comparative 
Example 2 


Reactor 


Novel Riser of 


Conventional 




the Present 


Riser 




Invention 




Density(20 °C), kgW 


745.3 


743.2 


Octane Number 






RON 


90.1 


89.8 


MON 


80.9 


79.8 


Induction Period, min 


800.0 


>485 


Existent Gum, mg/lOOml 


2.0 


2.0 


Sulfur, wt% 


0.01 


0.02 


Nitrogen, wt% 


0.003 


0.004 


Carbon, wt% 


86.51 


86.81 


Hydrogen, wt% 


13.42 


13.52 


Distillation, °C 






IBP 


48 


45 


10% 


75 


72 


30% 


97 


99 


50% 


118 


120 


70% 


144 


145 


90% 


173 


174 


EP 


203 


202 


Gasoline Composition, wt% 






Paraffins 


47.87 


14.41 


In which Iso-Paraffins 


43.86 


9.51 


Naphthenes 


7.45 


9.50 


Olefins 


20.51 


56.89 


Aromatics 


24.17 


19.20 



8. I have found out, upon scientific calculation, that if "DR2: DRi 
is equal to approximately 3: 1", the ratio of the velocity at the outlet of 
the first reaction zone to the velocity at the inlet of the second reaction 
zone, i.e. Ui/uq, will be 18.7. However, it is generally accepted in the art 
that the ratio of the velocity at the outlet of an iso-diameter riser( i.e. the 
so-called first reaction zone) to the velocity at the inlet of the 
iso-diameter riser( i.e. the so-called first reaction zone) is less than 3, 
and absolutely no more than 10. Thus, the Examiner's conclusion 
cannot be fulfilled. 

The technical background, prerequisites for deduction and detailed 
deduction procedure are described hereinbelow. 

Technical Background of the Patent Application 

The fluidized catalytic cracking process was developed to meet the 
high reactivity of a zeolite molecular sieve, and namely, the reaction 
time in the reactor is not supposed to be too long to produce undesirable 
products, such as excess coke and dry gas, and is generally controlled 
within 2 to 3 seconds. In order to meet the requirement of a short 
reaction time, isodiametric riser reactors disposed vertically are used for 
industrial purposes. According to the principles of fluidization, the 
effluent velocity at the inlet of the riser reactor is supposed to be 4 to 5 
m/s under normal operating conditions; otherwise, the device will not 
work. Cracking is a reaction resulting in an increased number of 
molecules. Such molecular expansion directly results in that the effluent 
velocity in the riser reactor increases sharply, and the effluent velocity 
at the outlet of the riser reaches as high as 15 to 18 m/s. Moreover, 
judging from the industrial design and application in the prior art, the 
diameter of the riser is not allowed to be too big. The diameter ratio of 
the upper and lower reaction zones in a riser is generally less than 1.5: 1 
in the art. 

During examination of this application, the Examiner in charge of 
this application drew the conclusion on the basis of the European patent 
EP0171460, after manually making a measurement in Figure 8 of 
EP0171460, that "DR2: DR] is equal to approximately 3: 1". In my point 
of view, this is impossible because of the following reasons: 
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Since reactions basically including cracking occur at the bottom of 
the riser (the so-called first reaction zone) after a hydrocarbon feedstock 
contacts with a hot regenerated catalyst, the conversion of the feedstock 
in said zone increases rapidly, and a large amount of small-molecular 
products are produced, which is the manifestation of molecular expansion. 
The reactions at the top of the riser (the so-called second reaction zone) 
basically include conversion reactions, such as hydrogen transfer and 
alkylation, and cracking does not occur there in a large scale. 
Consequently, the increase of the conversion of the feedstock slows down 
in said zone. These are all general knowledge in the field of catalytic 
cracking, as seen from, for example, Figure 1 which illustrates the 
distribution of conversion along the axial direction of the height of the 
riser (In-Su Han, Chang-Bock Chung, Chemical Engineering Science 56, 
2001, page 1954, see Annex I). 



1^' reaction zone 2"** reaction zone 




Riser height {m) 

Fig 1 : Distribution of conversion along the axial direction of the height of the riser 

It is generally accepted in the art that the ratio of the velocity at the 
outlet of an iso-diameter riser( i.e. the so-called first reaction zone) to the 
velocity at the inlet of the iso-diameter riser( i.e. the so-called first 
reaction zone) is less than 3, and absolutely no more than 10. This is 
exemplified by the following two examples: 

(1) Figure 2 illustrates the distribution of velocity along the axial 
direction of the height of the riser (In-Su Han, Chang-Bock Chung, 
Chemical Engineering Science 56, 2001, page 1978, see Annex II). It can 
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be seen that the ratio of the velocity at the outlet of the first reaction zone 
to the velocity at the inlet of the first reaction zone (i.e. Ui/uo in Figure 2) 
is less than 2. 



1^' reaction zone 2"'' reaction zone 




Riser height (m) 



Fig. 2 Distribution of velocity along the height of the riser 

(2) It can be readily seen from the reference Process Calculation 
and Technical Analysis of Catalytic Cracking (written by Cao Hanchang, 
Hao Xiren and Zhang Han, see Annex III, and the translation thereof, 
Annex IV) attached to this Declaration that, for example, the calculation 
results on page 242 (see also pages 5 and 6 of Annex IV) show a 
velocity of 11.45 m/s at the outlet of the riser, a velocity of 5.78 m/s at 
the inlet of the riser and a ratio of the former to the latter of 1 .98. 

In the Office Action issued by the Examiner ("DATE MAILED: 
01/23/2006"), the Examiner makes the following comments: 

Looking only to the data provided on page 49, lines 8-24, of Kmecak, 
et al (in the "specific embodiment" for FIG. VIII), the following can be 
said of the riser reactor: 

-riser reactor height^about 49 meters 

-average linear velocity of the suspension in the riser reactor=about 
24m/s. 

-total time required to traverse the riser reactor=about 2 seconds 
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-total time within the prelift zone of the riser reactor (i.e., the dry 
gas-steam-catalyst suspension residence time before contact with the 
atomized oil feed)= a fraction of a second up to 0.5 seconds 

-total time within the first reaction zone, second reaction zone, and 
first conjunct section, combined (i.e. a hydrocarbon residence contact 
time with catalyst particles) =up to about 1 or 1.5 seconds 

In the case that the total residence time within the riser reactor is set 
at 2 seconds, we can properly assume a residence time of roughly 0.5 
seconds within the prelift zone and a residence time of roughly 1.5 
seconds within the first reaction zone, second reaction zone, and first 
conjunct section. 

Looking now to Figure VIII, we can further assume from the 
proportions given in the drawing the riser reactor is roughly configured 
as follows : 

-Pre-lifi zone height (0. 1 m) =4. 9 m 

-First reaction zone height (0.3^H)=14. 7 m 
-second reaction zone height (0.5*H)=24.5 m 
-First conjunct section height (O.I *H)=4.9 m 

-the average linear velocity of the suspension in the prelift zone 
=9.8 m/s 

My deduction is based on the Examiner's analysis mentioned 
above and Figure 3 as shown hereinbelow: 



I Riser geometry | [ Local velocity | [ Average velocity | 




Uo=9.8m/s 



Fig. 3 
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The symbols are defined as follows: 

u — velocity, m/s; 

h — height of the riser, m 

t — residence time, s 

V — volume flow rate, m^/s; 

A — cross-sectional area of the riser, m^; 

d — diameter of the riser, m; 

X — ratio of the diameter of the second reaction zone to the diameter of 

the first reaction zone, dimensionless; 

p — pressure. Pa; 

M — molar mass, g/mol; 

R — gas constant, 8.314 J/mol»K; 

T — temperature, K 

p — density, kg/m^. 

Symbols t, p and V are defined by the following equations: 

t = !^ (1) 

U ^ ^ 

PxM 

v = uxA = uxTi:x (^"j (3) 

Prerequisites or Assumptions for said Deduction 

(1) No chemical reaction occurs in the prelift zone, the fluid is in a 
uniform motion state, and the average velocity is represented by uq. 

(2) The velocity in the first reaction zone exhibits a logarithmic 
increase tendency along the axial direction of the riser, and the average 
velocity in the first reaction zone is represented by the logarithmic mean 
value of the velocity at the inlet (uq) and the velocity at the outlet (ui), 
which is calculated as described hereinbelow: 

Ui - Up 

InUj — lnU(j 
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(3) The basic reactions in the second reaction zone are conversion 
reactions, and molecular expansion may be neglected. The state of the 
fluid may approximately be regarded as uniform motion state. The 
average velocity is presented by U2. 

(4) The velocity in the conjunct section between the first and 
second reaction zones exhibits a linear decrease tendency along the 
axial height of the riser, and the average velocity is presented by the 
mean value of the velocity at the inlet (ui) and the velocity at the outlet 
(U2), which is calculated as described hereinbelow: 

2 

Detailed Deduction Procedure : 

Based on the above-mentioned prerequisites or assumptions for 
deduction and the data disclosed in EP0171460 (the height of the prelift 
zone is 4.9 m, the height of the first reaction zone is 14.7 m, the height 
of the conjunct section between the first and second reaction zone is 4.9 
m, the height of the second reaction zone is 24.5 m, the total height of 
the riser is 49 m and the average velocity in the whole riser is 24.5 m/s), 
the residence time in the whole riser is calculated as described 
hereinbelow: 

4.9m 14.7m 4.9m 24.5m _ 49m 

~+ u^-p^-^ u^JHu,+-^ = 24.5m/s ^'^^ 
Inuj - Inuo 2 

It can be seen from the data provided by the Examiner that uo=9.8 
m/s, and thus Equation (4) may be written as: 

4.9 m 14.7m 4.9m 24.5m 

9.8m/s - 9.8m/s ^l±^i^ ^ 

InUi - ln9.8m/s 2 

According to the law of conservation of mass, the mass flowrates at all 
cross-sections of the riser are equal, and, in terms of Figure 3, the mass 
flowrate at the inlet of the conjunct section between the first and second 
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reaction zones (shown by the subscript "1") and the mass flowrate at the 
outlet of the conjunct section between the first and second reaction 
zones (shown by the subscript "2") are equal, as shown by the following 
equation: 

p, X V, = X (6) 



The general knowledge of the principles of catalytic cracking 
teaches us that the reactions occurring at cross-section "1" and 
cross-section "2" are basically conversion reactions and cracking does 
not occur in a large scale. Accordingly, it may be concluded that the 
number of molecules at said two cross-sections remains unchanged. 
Since the pressures and temperatures at said two cross-sections do not 
change much, it may be concluded, on the basis of the density definition, 
i.e. equation (2), that the densities at said two cross-sections are equal: 

= Pi (7) 

Equation (6) can thus be written as: 
V. = \\ (8) 



Upon substitution of the mass flowrate definition, i.e. Equation (3), 
into Equation (8), a new equation is obtained, which is as follows: 



U2 X n X (y) = Uj X n X (9) 



Equation (9) is transformed to obtain the following equation: 



The following monadic nonlinear equation about Ui is thus 
obtained by substituting Equation (10) into Equation (5): 

14.7m 4.9111 24.5m 

' — ar-^-Tn— = 2s (u) 



"i - 9.8m/s ' 
Inui - ln9.8m/s 



The Examiner believes that the ratio of the diameter of the second 
reaction zone to the diameter of the first reaction zone is 3, which 
means that x=3 in Equation (11). The following result is then obtained 
by solving Equation (11): 

Ui = 183.3m/s (13) 



Obviously, 

Uj _ 183.3m/s 
Uo ~ 9.8m/s 



= 18.7 (14) 



The ratio of 18.7 is far greater than the generally recognizable ratio 
of less than 3 in the art, and is even greater than 10. This is highly 
contrary to the fact pointed out in Item 8 herein that "it is generally 
accepted in the art that the ratio of the velocity at the outlet of an 
iso-diameter riser( i.e. the so-called first reaction zone) to the velocity at 
the inlet of the iso-diameter riser( i.e. the so-called first reaction zone) is 
less than 3, and absolutely no more than 10". 

Thus, the Examiner's opinion that "DR2: DRj is equal to 
approximately 3 : 1" is incorrect and impossible. 

9. I further declare that all statements made herein on personal 
knowledge are true and that all statements made on information and 
belief are believed to be true; and that these statements are made with the 
knowledge that willfully or intentionally false statements made by 
declaration to the United States Patent and Trademark Office are 
punishable by fine or imprisonment, or both, under 18 U.S.C. §1001; and 
that such willfully or intentionally false statements may jeopardize the 
validity and/or enforceability of the patent application with Appl. No. 
09/553,990 and any United States Patent issuing therefrom. 

Signed at D&^IV^^ this t tfAx day of 1[M>€ . 2007. 

Dr. Xu Youhao 
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Abstract 

The purpose of this study is to develop a detailed dynamic model of a typical FCC unit that consists of the reactor, regenerator, and 
catalyst transfer lines. A distributed parameter model is presented for the reactor riser to predict the distributions of the catalyst and 
gas-phase velocities, the molar concentrations of 4-lump species, and the temperatures. The regenerator is also modeled in such detail 
that the two-regime, two-phase behavior of typical fluidized beds can be accounted for. An efficient model solver was constructed on 
the basis of a modular approach in which the model equations are grouped into 12 modules each corresponding to a specific part of 
the unit and type of equations. The dynamic model is implemented in Part II of this paper for extensive simulation study on FCC 
processes. © 2001 Published by Elsevier Science Ltd. 

Keywords: Dynamic simulation; Modeling; Fluidized catalytic cracking; Fluidized bed; Modular approach 



1. Introduction 

A. fluidized catalytic cracking (FCC) process is a unit 
that converts heavy distillates like gas oils or residues to 
gasoHnes and rniddle distillates using cracking catalyst. 
Various types of FCC units are being operated over the 
world, including the Flexicracker and MODEL IV unit 
of Exxon, the Orthoflow unit of M.W. Kellogg, the FCC 
resid cracker of Stone & Webster. Since a typical FCC 
unit can process a large amount of the feedstock into 
more valuable products, the overall economic benefits of 
a refinery could be considerably increased if proper con- 
trol and optimization strategies are implemented. The 
first step toward the goal will be to develop an accurate 
model that represents the dynamics of the process. 

FCC processes are known to be very difficult to model 
because of the large process scale, complicated hydro- 
dynamics, and complex kinetics of both cracking and 
coke burning reactions. Accordingly, numerous papers 
on FCC modeling have had widely varying scopes and 
levels of modeling rigor. Some papers are focused on 
specific parts of the unit while others are only concerned 
with reaction kinetics or steady-state behavior. Several 
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studies on dynamic modeling of the whole FCC unit also 
have limited applicability due to the outdated process 
type or the oversimplified assumptions introduced in the 
modeling. 

Early papers on FCC modeling are mainly concerned 
with the reaction kinetics or steady-state process behav- 
ior. Weekman and coworkers (Jacob, Gross, Voltz, 
& Weekman, 1976; Voltz, Nace, & Weekman, 1971; 
Weekman & Nace, 1970) proposed the three-lump and 
10-lump cracking reaction models. Several variations of 
these models are also available in the Kterature. Weisz 
and Goodwin (1963, 1966a) and Weisz (1966b) carried 
out extensive studies on the kinetics of coke burning 
reactions in the regenerator. Hanc, Nakashio, and 
Kusunoki (1975) and Morley and de Lasa (1987, 1988) 
also estimated the kinetic parameters of the coke burning 
reactions for zeolite catalyst. De Lasa and Grace (1979), 
De Lasa, Errazu, Barreiro, and SoHoz (1981) and Errazu, 
de Lasa, and Sarti (1979) performed steady-state analyses 
of the regenerator using their models based on the two- 
phase theory. Elshishini and Elnashaie (1990a,b), Elshishini, 
Elnashaie, and Alzahrani (1992), Arbel, Huang, Rinard, 
Shinnar, and Sapre (1995a) and Arbel, Rinard, Shinnar, and 
Sapre (1995b) reported' the multiple steady states of an 
FCC imit found in their steady-state simulations. 

Dynamic models have been developed in several stud- 
ies for certain parts of the FCC process. Faltsi-Saravelou 
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and Vasalos (1991a) and Faltsi-Saravelou, Vasalos, and 
Dimogiorgas (1991b) developed a fiuidized-bed model 
for a regenerator comprising the dense bed and the free- 
board. The model is based on the two-phase theory 
(Kunii & Levenspiel, 1991) and contains some empirical 
equations obtained from fluidization experiments. 
Theologos and Markatos (1993) and Theologos, Nikou, 
Lygeros, and Markatos (1997) presented a quite detailed 
dynamic model of a reactor riser. Their model consists of 
a set of three-dimensional partial differential equations 
that describe the distributions of catalyst and gas vel- 
ocities, pressure, component concentrations, and temper- 
ature. The model appears to represent the maximum 
conceivable rigor for the riser dynamics, but it would be 
a formidable task to extend the rigor to other parts of an 
FCC unit as well. 

Dynamic models for the whole FCC unit have been 
presented in recent papers. Elnashaie and coworkers 
(Elnashaie & Elshishini, 1993; Elnashaie, Abasaeed, 
& Elshishini, 1995) extended their steady-state model to 
a simple dynamic model, and investigated the sensitivity 
and stability of a bed-cracking type FCC unit. Lopez- 
Isunza (1992) presented a distributed parameter dynamic 
model for mass and energy balance, but their model 
neglects the hydrodynamic aspects of the FCC unit. 
Moreover, their model did not take into account the 
mixing characteristics of catalysts in the regenerator. 
Mcfarlane, Reineman, Bartee, and Georgakis (1993) pre- 
sented quite a comprehensive model that covers most 
parts of a Model IV type FCC unit, including the reactor, 
regenerator, bowers, U-bends, compressors, furnace, and 
valves. But their models seem to be so oversimplified that 
the distributed characteristics of the reactor riser or the 
regenerator are not properly accounted for. Further- 
more, the numerous empirical equations in the model are 
tailored to the Model IV FCC unit, making it difficult to 
extend to other FCC units of different type. Arbel et al. 
(1995a,b) developed a model that can describe both the 
steady-state and dynamic behavior of an FCC unit being 
operated in both partial and full combustion modes. The 
model seems to be very successful to predict the steady- 
state behavior of an FCC unit such as steady-state and 
input multiplicities. Ali and Rohani (1997) presented 
a dynamic model that consists of several ordinary differ- 
ential equations. They presented analytic solutions of the 
differential equations after adopting pseudo steady-state 
assumptions. However, it is impossible to predict the 
after burning reaction because they do not take into 
consideration the freeboard region of the regenerator in 
their model. 

The purpose of this study is to develop a detailed 
dynamic model of a modern riser -type FCC unit that 
consists of the reactor, regenerator, and catalyst trans- 
port lines with slide valves. Once the dynamic model and 
a simulator implementing the model are available, they 
can be used as a tool (i.e., as a hypothetical plant) for 



carrying our various process system studies Hke control 
and optimization. The dynamic model is developed to 
complement the previous models presented in the litera- 
ture and is expected to serve as a tool for various process 
system studies on an FCC process. Hydrodynamic de- 
scriptions for the crucial parts of the unit are incorpor- 
ated into the model. Special attention has been paid to 
the reactor riser to predict the velocity distributions of 
the catalyst and gas phases, the molar concentrations of 
4-lump species, and the temperatures. The regenerator is 
also modeled in such detail that the two-regime (dense 
bed and freeboard), two-phase (emulsion and bubble) 
behavior of typical fluidized beds can be described. The 
cyclones in both the reactor and the regenerator are 
modeled to take into account their effect on overall 
dynamics of an FCC process. Pressure and valve equa- 
tions are included in the modeling of the catalyst trans- 
port Unes to describe of the interaction between the 
reactor and the regenerator. The model solver was con- 
structed on the basis of a modular approach. The model 
equations are grouped into 12 modules each of which 
corresponds to a specific part of the unit and type of 
equations, and then an efficient iterative scheme is em- 
ployed for convergence of all the modules. The results of 
steady-state and dynamic simulations using the de- 
veloped model are presented in Part II of this paper. 
Since the most important model fidelity consideration for 
FCC modeling is how well the model reproduces the 
steady-state behavior (Arbel et al., 1995a), we compared 
in Part II the steady-state behavior predicted by our 
model with those in the literature. 



2. Process modeling 

The FCC units can be grouped into a side-by-side type 
and a stacked type according to the configuration of the 
reactor and the regenerator. Fig. 1 shows a typical side- 
by-side type FCC unit that comprises a reactor and 
a regenerator joined by two catalyst transport lines. It 
accounts for large part of the modern-type FCC units 
being operated over the world. In the stacked type FCC 
unit the reactor is placed directly above the regenerator 
to straighten the catalyst transport line. Because the 
major physicochemical phenomena taking place in both 
the side-by-side type and the stacked type units are very 
similar to each other, we describe the modehng proced- 
ure only for the side-by-side type unit. 

Before presenting the detailed model of each process 
unit, we describe the major physicochemical phenomena 
taking place in the FCC unit. The cracking reaction of 
hydrocarbon feeds takes place in the riser, while the 
regeneration of deactivated catalyst by coke deposition 
takes place in the regenerator. The hydrocarbon feed like 
vacuum gas oil is injected into the bottom of the riser 
along with a small amount of steam (0.5-3 wt% of feed), 
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Fig. 1. Schematic diagram of a side-by-side type FCC unit. 

which leads to good atomization and reduces coke 
formation. The feed is subsequently vaporized upon con- 
tacting the hot catalyst flowing from the regenerator. The 
hydrocarbon vapors go through endothermic catalytic 
cracking reactions on their way up through the riser. 
Lighter hydrocarbons are produced as main cracking 
products along with by-product coke, which deposits on 
the catalyst surface and lowers the catalyst activity. The 
residence time of the catalyst and hydrocarbon vapors in 
the riser is typically in the range of 2-5 s. The riser top 
temperature is typically between 750 and 820 K and is 
usually controlled by regulating the flow of hot regen- 
erated catalyst to the riser. The disengaging section of 
a modern riser-type reactor only serves to separate cata- 
lyst particles from vapors. The product vapor from the 
disengaging section enters a main-fractionator where va- 
por products are separated into various boiling point 
fractions. The spent catalyst is separated from the vapor 
in the reactor cyclone and falls into the stripping section 
where the hydrocarbons remaining on the surface are 
removed by stripping steam. The stripped spent catalyst 
is recycled through a catalyst transport line to the regen- 
erator. 

In the regenerator which is operated in the fluidization 
regime, the coke is burnt ofl^the catalyst surface by the air 
blown into the bed. This combustion reaction serves to 
reactivate the catalyst and to maintain the bed temper- 
ature (950-980 K for a gas oil cracker, 980-1080 K for 
a resid cracker) high enough to supply the heat required 
for the vaporization and cracking reaction of the feed in 
the reactor. The regenerated catalyst flows continuously 
into the riser bottom through another catalyst transport 
line. The catalyst circulation rate between the reactor and 
the regenerator is controlled by the two slide valves 
attached on the catalyst transport lines. 



In this paper, the modeling is concentrated on the 
reactor, regenerator, and catalyst transport lines. The 
catalytic cracking unit has several auxiliary units: the air 
blower, expander, wet gas compressor, feed pre-heater, 
air heater, catalyst cooler, CO boiler, and so on. Because 
different catalytic cracking units have different designs 
for these auxiliary units, the modeling of these auxiliary 
units is not included in this paper. In the subsequent 
sections, we describe the structures and features of our 
model and modeling assumptions. The entire model 
equations are presented in detail in the Appendix. 

2.1. Reactor 

Fig. 2 shows a schematic of the reactor which is dissec- 
ted into several sections (the feed vaporization section, 
the riser, the disengaging-stripping section, and the reac- 
tor cyclones) for modeling purpose. The rectangles in the 
figure denote the control domains on which the model 
equations are derived, and the arrows denote flow direc- 
tions between sections. In our study the riser was 
modeled as a tubular reactor consisting of the catalyst 
and gas phases and the feed vaporization section as 
a junction of two streams. The disengaging-stripping 
section and cyclones are modeled as a continuous stirred 
tank (CST). Quasi-steady-state models are applied to the 
feed vaporization section, the riser, and the reactor cyc- 
lones because these sections or units (with the time con- 
stant in the range of a few seconds) show much faster 
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Fig. 2. Schematic diagram of the reactor model. 
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dynamics than other units like the disengaging-stripping 
section (with the time constant in the order of a few 
minutes) or the regenerator (with the time constant in the 
range of several minutes to an hour). Thus, the overall 
dynamics of the reactor can be reasonably taken to be 
dominated by the dynamics of the disengaging-stripping 
section. In fact, the quasi-steady state assumptions for the 
riser are common in the modeling of the modern-type 
FCC units (Elnashaie & Elshishini, 1993, 1995; Arbel 
et al, 1995a; Ali & Rohani, 1997; Ali, Rohani, & Corriou, 
1997). 

2.1.1. Feed vaporization section 

The feed vaporizes when it is mixed with the regen- 
erated catalyst in the feed vaporization section located at 
the bottom of the riser. The feed vaporization section is 
modeled as a pseudo-heat transfer system in which two 
streams (catalyst and feed) join. The temperature, pres- 
sure, and velocity of the vapor and the catalyst coming 
out of the vaporization section are calculated. These 
variables depend on the process variables such as feed 
temperature, feed characteristics, feed droplet size, cata- 
lyst temperature, and pressure. The volume expansion 
and temperature variation caused by the vaporization of 
liquid feed are considered in the modeling of the feed 
vaporization section. 

Because modern FCC units have a number of feed 
nozzles designed to help rapid vaporizations of the feed, 
it is possible to assume that the feed vaporizes instantly 
upon contacting the hot catalyst. This assumption is 
supported by the paper of Ali et al. (1997) in which it 
takes about 0.1 s (corresponding to about 3% of the 
residence time of the gas-catalyst mixture in the riser) to 
fully vaporize the feed. Based on the characteristics stated 
above, a quasi-steady-state model is postulated for the 
feed vaporization section. It is further assumed that the 
vapor feed is at the saturated state after vaporization and 
that the riser is operated under adiabatic condition. 

2.1.2. Reactor riser 

The riser is modeled as a one-dimensional tubular 
reactor without radial and axial dispersion. Momentum, 
mass, and energy balance equations are derived for both 
the catalyst phase and the gas phase under the following 
assumptions: (1) the cracking reactions only take place in 
the riser; (2) the riser wall is adiabatic; (3) dispersion and 
adsorption inside the catalyst particles are negligible; (4) 
the coke deposited on the catalyst does not affect the 
fluid flow; (5) the coke has the same properties as the 
catalyst; (6) the riser dynamics is fast enough to justify 
a quasi-steady-state model. 

The momentum equations are included in the riser 
modeling to describe the significant variation of the lin- 
ear velocities of the catalyst and gas phases along the 
riser. There is significant molar expansion with increas- 
ing conversion of the feed and consequent acceleration of 
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Fig. 3. Profiles of gas oil conversion calculated with and without 
momentum equations. 
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Fig. 4. Four-lump model for gas oil cracking reactions. 



both the catalyst and gas phases (Theologos & Mar- 
katos, 1993; Theologos et al., 1997). But, most studies on 
the riser modeling in the literature assumed the linear 
velocity to be constant throughout the riser, probably to 
keep the momentum equations from making the entire 
model stiff and complex. To investigate the degree of 
error arising from assuming constant velocity, we com- 
pared the gas oil conversions obtained by including and 
excluding the momentum equations in the riser model, 
respectively. Fig. 3 shows that the conversions at the riser 
exit can deviate from each other by as large as 9.6%. It 
becomes clear that the momentum equations can explain 
the molar expansion and catalyst slip in the riser and can 
considerably reduce the error arising from assuming con- 
stant values for such properties of the catalyst and gas 
phases as velocities, volume fractions, density, and pres- 
sure. Furthermore, the momentum equations enable one 
to predict the pressure drop across the riser and conse- 
quently the pressure in the feed vaporization section. 

The mass balance equations are used for predicting the 
component weight fractions of the four lumps: gas oil, 
gasoline, light gases (C1-C4), and coke. Fig. 4 shows the 
four-lump model of Lee, Chen, and Huang (1989) that 
expresses the kinetics of cracking reactions in the riser. 
This model which is a revised version of the three- 
lump model of Weekman and Nace (1970) takes the coke 
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yield into account. The more detailed 5-lump model 
(Ancheyta-Juarez, Lopez-Isunza, Aguilar-Rodriguez, & 
Moreno-Mayorga, 1997) or 10-lump model (Jacob et al., 
1976) may be used if one wants to predict the yields of 
subdivided fractions. It is assumed that cracking of gas 
oil is a second-order reaction but that of gasoline is 
a first-order reaction, and that the reactions take place 
only in the gas phase. Since the overall cracking rate is 
affected by the catalyst activity, its effect should be in- 
corporated into the kinetic expressions. However, the 
catalyst activity is very hard to model with certainty since 
the intrinsic catalyst activity and the coking rate depend 
upon feed composition, the addition of catalyst activity 
modifiers, and the rate at which the catalyst is replaced 
with fresh catalyst. In this study it is assumed that coke 
deposition on catalyst is the major cause of catalyst 
deactivation, and that its effect can be described by 
a catalyst deactivation function of the following form: 

where the catalyst deactivation coefficient is related to 
the temperature and the feedstock composition by the 
following equation: 

Heavier gas oil feeds generally contain large amount of 
Conradson carbon and may affect the catalytic cracking 
rate. It is assumed that the Conradson carbon in the 
virgin feedstock is directly converted to coke at the en- 
trance of the riser and deactivates the catalyst. 

The catalyst temperature is higher than the gas tem- 
perature throughout the riser and the consequent heat 
transfer provides the heat required for the vaporization 
and endothermic reaction of the feedstock. The heat 
transfer between the two phases and the heat of cracking 
reaction are considered in the energy balance for the riser 
to predict the temperature of each phase. 

2.1.3. Disengaging-stripp ing section 

The disengager and the stripper of the reactor are 
combined into a single section called the "disengaging- 
stripping section" in our model. The section is modeled 
as a perfectly mixed continuous tank with no reaction 
taking place. Because the catalyst is immediately separ- 
ated from the product vapor via the cyclones, further 
cracking reaction seldom occurs in the disengaging-strip- 
ping section. The disengaging-stripping model compris- 
ing the coke, catalyst, gas component, energy, and 
pressure balances is used to calculate several major state 
variables: the coke on catalyst after stripping, the catalyst 
and gas holdups, the concentrations of the gaseous com- 
ponent into the main-fractionator, and the reactor tem- 
perature and pressures. 



A portion of hydrocarbons is trapped in the pores of 
the catalyst in the riser. The unremoved hydrocarbons 
are converted in the regenerator to coke (called catalyst- 
to-oil coke or occluded coke) and thus cause the loss of 
product yield and an increase in the air required for 
catalyst generation. Some of the hydrocarbons can be 
removed in the stripper by steam stripping the rate of 
which depends on the catalyst-to-oil ratio and steam 
injection rate. In this study, the following exponential 
type stripping function is developed to account for the 
catalyst-to-oil coke: 

CckST = CcisTO + Kso expj^ - .^iffsM^k J. (3) 

In the modeling of an FCC unit, pressure should be 
treated as one of the most important process variables 
because the pressure difference between the reactor and 
the regenerator determines the catalyst circulation rate 
and thus the degree of interaction between the two 
vessels. The pressure in the disengaging-stripping section 
is calculated from the ideal gas law and the pressure at 
the bottom of the reactor is estimated using the static 
head exerted by the stacked catalyst particles. 

2.1.4. Reactor cyclones 

Both the reactor and the regenerator of an FCC unit 
are usually equipped with several multi-stage cyclones to 
separate catalyst particles from entraining vapors. All the 
cyclones in the reactor are lumped into one modeling 
unit, which is then described as a continuous stirred tank 
(CST). It is assumed that aU the reactor cyclones have the 
same dimensions and that the catalyst holdup in each 
cyclone is equal to each other. It is further assumed that 
the flow rate of the catalyst exiting the cyclone dip-legs is 
only determined by the catalyst holdup in the cyclones. 
Because no cracking reaction takes place in the reactor 
cyclones, the cyclones are in thermal equilibrium with the 
disengaging-stripping section. Therefore, only the mass 
balances are required to calculate the following state 
variables: the reactor cyclone inlet velocity, the catalyst 
holdup in the reactor cyclones, and the mass flow rate of 
the fluids exiting the reactor. 

2.2. Regenerator 

The regenerators of most FCC units are operated as 
a fluidized bed which consists of two sections commonly 
called a dense bed and a freeboard. The dense bed con- 
tains the majority of the catalyst in the regenerator and 
thus accounts for the generation of spent catalyst by coke 
burning reactions. The well-known two-phase theory of 
fluidization (Kunii & Levenspiel, 1991) describes the 
dense bed as comprising a particle-lean bubble phase and 
a particle-rich emulsion phase surrounding the bubbles. 
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Fig. 5. Schematic diagiam of the regenerator model. 



The gas bubbles erupt at the surface of the dense bed and 
eject the catalyst particles into the freeboard. Since the 
catalyst bulk density is very low in the freeboard, the 
coke burning reaction is slower than in the dense bed, but 
there take place the after-burning reactions that oxidize 
carbon monoxide to carbon dioxide. The catalyst par- 
ticles entrained in the stack gas are separated in the 
cyclone and then are recycled to the dense bed. The 
physical situations described here are illustrated in 
a schematic diagram of the regenerator shown in Fig. 5. 
The dense bed is modeled as a hybrid reactor that em- 
ploys a mixed-tank model for energy and coke balances 
but a tubular reactor model for gas component balances. 
The freeboard is modeled as a two-phase tubular reactor, 
and the cyclone as a mixed-tank. 

Because the regenerator contains much larger amount 
of catalyst than the reactor, the residence time of catalyst 
in the former (usually falling within the range of several 
minutes to an hour depending on the type of FCC units) 
is much longer than in the latter. Thus, the overall dy- 
namics of an FCC process is largely dominated by the 
dynamics of the regenerator, which necessitates quite 
a detailed description of the regenerator for reliable pre- 
dicting of the entire process behavior. In this study, 
dynamic balances are set up for the coke, catalyst, gas 
component, and energy in every dissected region. 

We begin with the reaction kinetics and the overall 
hydraulics of the regenerator. 



2.2.1. Coke burning kinetics 

The typical coke consists of various chemical com- 
pounds including hydrocarbons, sulfuric and nitro- 
genous compoimds. Since the hydrocarbons are 
dominant, however, in this study the coke is assumed to 
have the formula CH, and to burn according to the 
following reaction schemes (Weisz & Goodwin, 1963, 
1966a; Weisz, 1966b; Hano et al., 1975): 

CH, -I- (0.5 + 0.25^)02 ^ CO -I- aSgH^O 

I'll — l<^tRoCckiCo^i, 
CH, + (1 + 0.25q)O2 -+ CO2 + O.S^HzO 

r2, = k^RoC^^iCo^i, i = E,F}, (4) 
where the rate constants can be represented by 



1 - 



1 - 



(5) 



Besides the coke burning reaction, there takes place 
the oxidization reaction of the carbon monoxide. It may 
take the form of either homogeneous oxidation on the 
gas phase with the kinetics proposed by Howard, Willi- 
ams, and Fine (1973) 

2CO + O2-* 2CO2 rsi = k:,RoCcoiCZ^C?i^o■„ 

i = B,E,F (6) 

or heterogeneous oxidation in the presence of the catalyst 
which contains small amounts of oxidation promoters 
like vanadium, nickel, or copper. Krishna and Parkin 
(1985), and Ali and Rohani (1997) reported the following 
kinetics of heterogeneous oxidation: 



2CO H 



^2C02 



(7) 



2.2.2. Overall regenerator 

The total catalyst and gas holdups, superficial gas 
velocity, average gas density, and pressures are calculated 
from the balances around the regenerator as a whole. 
These values are then used in the model equations for 
each section. The dynamics of the total catalyst and gas 
holdups are important in describing the interactions be- 
tween the reactor and the regenerator because these 
variables strongly affect the catalyst flow rates between 
the two vessels. The superficial gas velocity varies along 
the regenerator due to the volume change. Since the mass 
influent to the gas phase increases in proportion to the 
coke burned off the catalyst, the superficial gas velocity 
can be correlated to the average gas density in the regen- 
erator: 



PgRO^Rl 



(8) 
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In the above equation, 5gRG denotes the relative extent to 
which the gas phase mass flow rate increases by the coke 
burning reaction. It is defined to vary from 0 at the 
bottom to 1 at the top of the regenerator, and can be 
calculated by the following expression: 



\j).5qpgRG 



-(0) ' 
HjOD J 



= D, F. 



(9) 



The pressures at the dense bed outlet and the bottom of 
the regenerator are calculated by applying average gas 
density and are then used to compute the flow rates 
through the catalyst transport lines. 

2.2.3. Dense bed 

In this study, the dense bed is modeled to consist of the 
emulsion and bubble phases on the basis of the two- 
phase theory. This contrasts with the approaches of 
Arandes and de Lasa (1992) and Arbel et al. (1995a) who 
modeled the bed to consist of gas and catalyst phases. 
The two-phase model can describe the whole range of 
fluidization regimes that cover the operating conditions 
of various types of regenerators. Our dynamic balances 
for the dense bed are based on the hybrid reactor model 
combining a mixed-tank reactor for coke and energy 
balances and a tubular reactor for gas component bal- 
ances. This also contrasts with Faltsi-Saravelou and 
Vasalos (1991), Ali and Rohani (1997), and AH et al. 
(1997) who used mixed-tank reactors for gas component 
balances in the emulsion phase. 

The following assumptions are introduced to facilitate 
the modeling of the dense bed: (1) the bubbles do not 
contain any catalyst particles; (2) all the catalyst particles 
ejected into the freeboard are returned to the dense bed 
via cyclones; (3) the catalyst particles in the dense bed are 
perfectly mixed due to the recycle from the cyclone as 
well as their own circulation in the bed; (4) the catalyst in 
the dense bed is in thermal equilibrium with the gases; (5) 
the gas densities of both the bubble and emulsion phases 
are the same as each other; (6) the gases in both the 
bubble and emulsion phases are in tubular flow and 
(7) the gas in the regenerator consists of oxygen, carbon 
monoxide, carbon dioxide, water vapor, and nitrogen. 

The hydraulic properties like the bed heights, phase 
volume fractions and velocities represent the major 
modeling elements for the regenerator because they are 
closely connected to the mass and energy balances in the 
bed. But the detailed description of the hydraulics is quite 
involved due to the inherent complexity of the fluidiz- 
ation phenomena. In this study, an empirical correlation 
between the volume fraction and the average superficial 
gas velocity was developed by a nonlinear regression of 
the experimental data of Schnitzlein (1987). 



Then, the bed height can be calculated on the basis of the 
above equation. The volume fraction of the bubble gas is 
predicted by the following expression proposed for fast 
bubbles and Geldart A type particles (corresponding to 
most FCC catalysts) (Davidson & Harrison, 1963; Kunii 
& Levenspiel, 1991): 



UgRG - 



(11) 



where the minimum fluidization velocity (u^^) and bubble 
rising velocity {vgg) can be estimated by empirical expres- 
sions in the literature (Kunii & Levenspiel, 1991). 

From the gas component, coke, and energy balances, 
the following state variables are computed: the molar 
concentrations of gaseous substances along the axial 
position of the dense bed, the coke on catalyst, and the 
dense bed temperature. In addition, the volume fractions 
and the dense bed height are calculated from other mass 
balances and empirical expressions. 

2.2.4. Freeboard 

The freeboard is modeled as a tubular reactor. Em- 
ploying a freeboard model is important in predicting the 
after-burning reactions in the freeboard. It is assumed 
that the catalyst in the freeboard is in thermal equilib- 
rium with the gaseous phase and is in plug flow. In most 
models for the freeboard in the literature, the hydraulic 
properties like the volume fraction of catalyst is usually 
assumed to be constant along the freeboard. But, the 
concentration of the catalyst tends to decay exponenti- 
ally with increasing height in the freeboard (Kunii 
& Levenspiel, 1991). In this study, the catalyst volume 
fraction in the freeboard is estimated with the following 
empirical equations (Li & Kwauk, 1980): 



I- {^cD - £?F)exp[ - p,{z - Hd)1 



(12) 



where e?f denotes the saturated carrying capacity of gas. 
We carried out nonlinear regression using the experi- 
mental data ofZenz and Weil (1958) to obtain the follow- 
ing correlation equation for e^p: 



_ PgRG 



(13) 



= 0.3418exp[ - 0.975 l£i^j5cj] -t- 0.1592. 



(10) 



The freeboard model comprising the gas component, 
energy, and coke balances is used to calculate several 
state variables: the molar concentrations of gaseous sub- 
stances, coke on catalyst, and the temperature along the 
axial position of the freeboard. 

2.2.5. Regenerator cyclones 

Since the catalyst particles collected by the cyclones 
are recycled to the dense bed, there is significant interac- 
tion between the dynamics of the two sections. The 
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regenerator cyclone is also modeled as a continuous 
stirred tank like the reactor cyclone. A pseudo-steady 
state is postulated because the catalyst residence time in 
the regenerator cyclone is much shorter than in other 
regenerator sections. The after-burning of CO to COj in 
the regenerator cyclones (including a plenum chamber) is 
likely to occur when there is an excess of both oxygen and 
carbon monoxide. In this study, the temperature rise 
across the regenerator cyclones is estimated using the 
following simple correlation (Kurihara, 1967; Hovd 
& Skogestad, 1993): 

TcY2 = + cr/Sr^'. (14) 

The regenerator cyclone model yields the mass flow 
rate of the catalyst into the regenerator cyclones, the 
catalyst holdup in the regenerator cyclones, and the re- 
generator cyclone temperature. 

2.3. Catalyst transport lines 

The catalyst circulation rate between the reactor and 
the regenerator is controlled by the two slide valves 
attached on the catalyst transport lines. The catalyst 
circulation rate through each catalyst transport line is 



Table 1 



Types of the model equations derived for each section 


Process units 


PDEs 


ODEs 


AEs 


No. of Eqs. 


Feed vaporization section 


0 


0 


7 


7 


Riser 


0 


9 


27 


36 


Disengagjng-stripping section 


I 0 


7 


14 


21 


Reactor cyclone 


0 


0 


5 


5 


Overall regenerator 


0 


2 


17 


19 


Dense bed 


12 


2 


44 


58 


Freeboard 


7 


0 


23 


30 


Regenerator cyclone 


0 


0 


5 


5 


Catalyst transport hnes 


0 


0 


8 


8 


Total No. of Eqs. 


19 


20 


150 


189 



determined by the pressure drop across a slide valve and 
can be calculated using the following simple valve equa- 
tion: 

f CSV = /cX(x„)APgif . (15) 

Because the catalyst residence time in the catalyst trans- 
port lines is in the order of a few seconds for stacked type 
FCC units and ten seconds for the side-by-side type FCC 
units, it is possible to ignore the transport lag throughout 
the catalyst transport lines in their model. The following 
variables are calculated using the catalyst transport line 
equations: the mass flow rates through the spent and 
regenerated catalyst transport lines, the coke on catalyst, 
and the catalyst temperatures. 



3. Numerical algorithms 

The entire dynamic model equations for the FCC 
process represent a mixed system of diff"erential (ordinary 
and partial) and algebraic equations. The total of 189 
equations comprising the model are arranged in the 
Appendix and classified in Table 1, The remainder of this 
paper describes our numerical algorithms for solving the 
model equations. The estimation methods for the para- 
meters in the model and the results of simulation runs are 
presented in Part II of this paper. 

3. 1. Degree of freedom analysis 

A degree of freedom analysis should be the first step in 
the development of a dynamic simulator, especially for 
a large and complex system like an FCC process. 
The FCC model presented in this study contains 202 
unknown variables as summarized in Table 2 when ex- 
cluding such determinable parameters as equipment di- 
mensions, physical properties, feedstock characteristics. 
Hence, the dynamic model is found to have 1 3 degrees of 
freedom ( = 202 variables — 1 89 equations), which 



Table 2 

Unknown variables in the model equations for the FCC unit 



Variable type 


Symbols 


No. of unknown variables 


Temperature 


T 


14 


Pressure 


P 




Holdup 


w, V. H 


12 


Density & concentration 


C,C,„f,Y, y,p 


51 


Volume fraction 




8 


Flow rate & velocity 


F,v,u 


30 


Heat rate 




18 


Reaction coefficient 


fc, R, r, a„ AH, (f>. 


43 


Valve stem position 




3 


Others 


Ci, Cf, dt,f, G, M„, Re, S, 


15 


Total No. of unknown variables 




202 



I.-S. Han, C.-B. Chung / Chemical Engineering Science 56 (2001) 1951-1971 



means that we have to provide 13 additional relations 
among the variables before trying to obtain a unique 
solution. Typically, these relations are provided by the 
control laws which connect the manipulated variables to 
the controlled variables under a given control configura- 
tion for the process. Since this paper is focused at present 
on the dynamic behavior of the FCC process in an 
open-loop mode, however, the 13 degrees of freedom are 
fixed by specifying the same number of variables appro- 
priately. The variables chosen for this purpose are 
presented in Part II of this paper where our process 
simulation results are presented. 

3.2. Model equation solvers 

The model solver was developed on the basis of 
a modular approach (Ponton, 1983; Hillestad & Hertz- 
berg, 1986) in which the equations were grouped into the 
modules and then all the modules were sequentially 
solved at the same time step. In this study the model 
equations are classified into modules, each of which cor- 
responds to a specific physical section of the reactor or 
regenerator. Then some modules are further divided 
according to the type of equations or the domain of 
analysis. Fig. 6 shows the overall structure of the model 
comprising twelve modules and the flows of information 
between them. As evident from the figure, the modules 
are severely coupled to each other due to the recycles and 
interactions between the sections. 

Each individual module was solved using a solver 
specifically chosen for the type of equations in the mod- 
ule. There are four types of equations, and each type was 
solved by the following numerical methods: 




Fig. 6. Information flow between the modules. 



(1) Nonlinear algebraic equations: The HYBRDl rou- 
tine of MINPACK (Garbow, Hillstrom, &. More, 1980) 
was used. It adopts the PoweU's hybrid algorithm based 
on a Newton's method. 

(2) Partial differential equations: The method of lines 
(MOL) was employed that converts the partial differen- 
tial equation into a set of ordinary differential equations 
using finite difference approximation of the spatial deriv- 
atives and then integrates the set of equations with re- 
spect to time. First-order backward difference was used 
to approximate the convection terms and second-order 
central difference to approximate the diffusion terms. 
A non-uniform grid scheme was also adopted to enhance 
the flexibility of meshing. 

(3) Coupled differential and algebraic equations: The 
Pezold-Gear's BDF algorithm (Petzold, 1982) was used 
to solve the momentum equation of the reactor riser that 
belongs to this class of equations. The algorithm 





Fig. 7. Modular iteration scheme for convergence. 
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uses backward difference formulas and shows powerful 
performance for highly stiff problems. 

(4) Ordinary differential equations: An LSODE routine 
(Hindmarsh, 1983) based on the GEAR'S algorithm was 
used. It implements a predictor-corrector method for 
non-stiff problems and a backward difference formula 
method for stiff problems. The routine was modified in 
this study to allow the concurrent integration of multiple 
modules. It also employs banded matrix algorithms to 
save storage and computation time. 

Fig. 7 shows the iterative procedure employed in this 
study to obtain convergent solutions for all the modules 
at every simulation time step. In the figure, n denotes the 
time step counter, and k the iteration counter at each 
time step, and x"'^' the vector of the dependent variables 
of the ith module at the kth iteration step and time n At. 
Each module is solved in turn at every iteration using an 
appropriate equation solver discussed previously. An im- 
pUcit solution in time is sought for in the manner analog- 
ous to the Gauss-Seidel algorithm for solving hnear 
algebraic equations. More specifically, in calculating the 
time-implicit vector x}"tVi'^', the updated iteration vec- 
tors Xi"fc+i^', ... ,x-"_'*ij^+'i of the preceding modules and 
the current iteration vectors xlVi!*"^') • • • > of the 

subsequent modules are treated as fixed parameters. On 
the other hand, when the ith module includes differential 
equations, the vector xf^' obtained at the previous time 
step is always used as the initial conditions for time 
integration. This cycle of computation is repeated 
until all the modules satisfy their respective specified 
tolerance, and then the simulation proceeds to the next 
time step. 



4. Conclusions 

A detailed dynamic model of a modern riser-type FCC 
unit that consists of the reactor, regenerator, and catalyst 
transport lines with slide valves was developed on the 
basis of the conservation principles. For modeling pur- 
poses the reactor was divided into the feed vaporization 
section, riser, disengaging-stripping section, and reactor 
cyclones, while the regenerator into the dense bed, free- 
board, and regenerator cyclones. Distributed parameter 
models were presented for the reactor riser, the dense 
bed, and freeboard of the regenerator, while lumped 
parameter models were postulated for the other parts of 
the unit. An efficient model solver was constructed on the 
basis of a modular approach in which the model equa- 
tions are grouped into 12 modules each corresponding to 
a specific part of the unit and type of equations. The 
dynamic model is implemented by a Fortran code and 
the simulation results are presented in Part II of this 
paper. 



Notarion 

All notations should be interpreted according to the 
following rules. The general form of a symbol is A^f/*, 
where N denotes the property of interest and consists of 
a letter with or without subscripts. The subscript i in 
lower-case represents substance or species and may have 
one or more letters. The subscript J in upper-case desig- 
nates the part or section of the FCC unit in which the 
substance i exists and also have one or more letters. The 
subscript k is optionally used to provide additional in- 
formation when necessary. Hence, for example, CpgRs de- 
notes the heat capacity of the gas in the riser. It should be 
noted that subscripts are often omitted for commonly 
used symbols with obvious meaning, e.g., which de- 
notes the catalyst density in all the sections of the FCC 
unit. 



A surface area, m^ 

effective heat transfer area between the 
dense bed and the surroundings, m^ 

Aig, Big, C,g cocfficients of the Antoine equation for 
gas oil feedstock 

Apjise effective interface heat transfer area per 

unit volume between the catalyst and gas 
phases in the riser, m^/m^ 

A^Te effective heat transfer area between the 

reactor and the surroundings, m^ 

Cj- temperature rising factor in the regener- 

ator cyclone equation, K 

C mole concentration, kg mole/m^ 

coke on catalyst, kg coke/kg catalyst 

CckSTo minimum coke content attainable by 

stripping, kg coke/kg catalyst 

Ca drag coefficient 

Cf catalyst-gas friction coefficient, kg/m^ s 

Cp heat capacity, kJ/kg K 

Cp mean heat capacity, kJ/kg K 

d average diameter, m 

rfcm minimum particle diameter that should 

be completely separated by a cyclone, m 

D diameter, m 

E activation energy, kJ/kg mole 

E^^ exponent in the stripping function, s/kg 

F mass flow rate, kg/s 

/ mole fraction or friction coefficient 

mole fraction at the bed inlet 

/„(x„) valve flow characteristic function: 

Xv/Lo'-v + (1 — a-j^u]""^ for a linear valve, 
and xj/[a„ + (1 - a„)x^]*'-* for an in- 
creasing sensitivity valve 

Fa mass flow rate of the air entering the 

dense bed, kg/s 

Fas mass flow rate of the dispersion steam 

entering the feed injection section, kg/s 
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1961 



G 



H,h 
k 



A^CYl 



mass flow rate of the liquid feed entering 
the feed injection section, kg/s 
mass flow rate of the stripping steam en- 
tering the disengaging-stripping section, 
kg/s 

acceleration due to gravity, 9.8 m/s^ 

stress modulus, kg/m s^ 

inlet height of each cyclone in the reactor, 

m 

inlet height of each cyclone in the regen- 
erator, m 

interface heat transfer coefficient between 
the catalyst and gas phases in the riser, 
kJ/m^ s K 

height of regenerator, m 
height of riser, m 
height of disengaging section, m 
height of stripping section, m 
specific enthalpy, kJ/kg 
heights, m 

rate coefficient of a reaction 
rate coeflficient of the coke burning reac- 
tion, m'/kg mole s 

frequency factor in the Arrhenius expres- 
sion for the coke burning reaction, 
m^/kg mole s 

rate coefficient of the homogeneous CO 
oxidation, m^/kg mole s 
frequency factor in the Arrhenius expres- 
sion for the homogeneous CO oxidation, 
m^/kg mole s 

rate coefficient of the catalytic CO oxida- 
tion, m^^'^/kg catalyst kg mole°'^ s 
frequency factor in the Arrhenius expres- 
sion for the catalytic CO oxidation, 
m^^'^/kg catalyst kg mole^"^ s 
rate coefficient of the four-lump cracking 
reaction, 1/s 

frequency factor in the Arrhenius expres- 
sion for the four-lump cracking reaction, 
1/s 

frequency factor in the stripping function, 
kg coke/kg catalyst 
valve flow rating factor, kg/s kPa°"* 
flow rating factors of the cyclone dip-legs, 
kg/s m^-^ 

interchange coefficient between the 
bubble and emulsion phases, 1/s 
length, m 

lengths of the catalyst transport lines, m 
molecular weight 

numbers of cyclones in the reactor and 

regenerator, respectively 

number of turns of 360° made spirally 

downward in a cyclone 

pressure, kPa 



^conv 



pressure at the main fractionator, kPa 
hydrocarbon to carbon atomic ratio in 
coke 

rate of heat generation or heat removal 
by reaction, kJ/s 

rate of heat loss to the surroundings, kJ/s 
reaction rate 

ideal gas constant, 8.3143 kPam^- 
kg mole K or kJ/kg mole K) 
aromatics-to-naphthenes ratio in bquid 
feedstock 
Reynolds number 

average sphericity of catalyst particles 
total mass interchange rate between the 
emulsion and bubble phases, 1/s 
temperature, K 

temperature of the air entering the dense 
bed, K 

temperature of the dispersion steam en- 
tering the feed vaporization section, K 
temperature of the surroundings, K 
temperature of the liquid feedstock enter- 
ing the feed vaporization section, K 
reference temperature, 298 K 
temperature of the stripping steam enter- 
ing the disengaging-stripping section, K 
superficial velocity, m/s 
average superficial velocity, m/s 
overall heat transfer coefficient between 
a process unit and the surroundings, 
kJ/m^ s K 

interstitial velocity, m/s 

superficial gas velocity at minimum 

fluidizing condition, m/s 

volume, m^ 

holdup, kg 

axial coordinate in the riser, m 

axial coordinate in the spent catalyst 

transport line, m 

axial coordinate in the regenerated cata- 
lyst transport line, m 
valve stem position, [0—1] 
weight fraction, [0 — 1] 
weight fraction of Conradson carbon 
residue in a feedstock 
gas oil conversion, [0—1] 
axial coordinate in the regenerator, m 
gas compressibility factor 



Greek letters 



catalyst deactivation coefficient 

preexponential factor of 

exponent for representing 

gas-phase heat diffusion coefficient, m^/s 
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valve head differential at maximum 


CYl 


reactor cyclones 




flow/valve head differential at zero flow, 


CY2 


regenerator cyclones 




[0 - 1] 


D 


dense bed 




catalyst decay constant, 1/s 


d 


cyclone dip-leg 


"a 


relative extent to which the g3^s phase 


ds 


dispersion steam 




mass flow rate increase by the coke burn- 


E 


emulsion phase 




ing reaction 




effective 




heat of formation, kJ/kg mole 


F 


freeboard (dilute phase) region 






FS 


feed vaporization section 




ing reaction modelj kJ/kg 








heat of stripping, kJ/kg 


% 


gasoline 




heat of vaporization of liquid feedstock in 




gas oil 




the feed vaporization section, kJ/kg 




light (C1-C4) gases 






H2O 






void fraction in the dense bed at min- 


hg 


hydrocarbon gases 










'/ CY 1 J '/ CY 2 


collection efficiencies of cyclones, fO — 1] 


Ig 


liquid feedstock (gas oil) 


CLl, CL2 


angles of inclination of the catalyst trans- 




maximum or minimum 




port lines, degrees 


MF 


main-fractionator 




viscosity, kg/m s 


N2 


nitrogen 




stoichiometry coefficient for the coke 


O2 


oxygen 




burning reaction 


RG 


regenerator 


n 




RS 






density, kg/m 


RT 


disengaging-stripping section (or reactor) 




intrinsic CO 2 /CO molar ratio in coke 








catalyst deactivation function 


ss 


stripping steam 


X 




SG 


stack gas valve 


Q 


cross sectional area m 


ST 


stripping section 




cross sectional area of each reactor cyc- 


SVl 


slide valve of the regenerated catalyst 




lone dip-leg, m^ 




transport line 


i^CYli 


cross sectional area of each reactor cyc- 


SV2 


slide valve of the spent catalyst transport 




cross^sectional area of each regenerator 






CY2d 


cross sec lona r ^ g 


wv 


water vapor 






wg 


wet gases 


"CY2i 


cross sectional area ol each regenerator 








^^'^ss" ecVc^^krarea of re enerator m^ 








cross sec lona area o regenera or, m 






^RS 


cross sec lona area o nser, m 


uperscrip s 










fl 

incoming ow 






out 




u scrip s 






atariaWe evaluated at v 






(0) 


a vTr^bleTTOhTafed^t'^the entrance of 








unit 




b^'^bbl 
u e p ase 




evaluated at the dense bed exit 




bottom 




evaluated at the top of the regenerator 




single bubble or catalyst bulk 




evaluated at the riser exit 


c 


catalyst 


it) 


evaluated at time f 


cc 


Conradson carbon 


(Td) 


evaluated at the temperature of the dense 


cfc 


coke 




bed 


CLl 


regenerated catalyst transport hne (to re- 


(Tf) 


evaluated at the temperature of the free- 




actor) 




board 


CL2 


spent catalyst transport line (to regener- 


(X) 


evaluated through the axial position of 




ator) 




the riser 


CO 


carbon monoxide 




evaluated through the axial position of 


CO2 


carbon dioxide 




the regenerator 
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, ^RS ^^cRS _ CfRs{VgRS — ^cRs)^RS 



(A. 8) 



Appendix A. Process model equations 

A.l. Feed vaporization section 

The feed vaporization section is modeled as pseudo- 
heat transfer system in which two streams (catalyst and 
feed) join. The catalyst temperature after feed vaporiza- 
tion is calculated by the following energy balance assum- 
ing adiabatic operation, and the vapor temperature is 
calculated by the Antoine equation. 



--^ ^^—LCpig{Tgfs — Tig) 

i'cCLl ^DC 



{T,ps - T,,) + AH, 



Alg - l0g(PFSy9 oPS) 



fRT + A^s = i' 



- ip'rI - n¥'\ 



(A.1) 
(A.2) 



(A.3) 
{A.4) 



The velocities of the vapor and the entrained catalyst 
entering the riser are calculated by the continuity equa- 
tions: 



P9FS(1 - 



Pc^cCLl^Rs' 



A.2. Reactor riser 



(A.5) 
(A.6) 



The riser is modeled as a one-dimensional tubular 
reactor without radial and axial dispersion. The mo- 
mentum equations of both the catalyst and gas phases in 
the riser are included in the riser modeling and can be 



dPgJiS ^ CIrs ^Frs _ CfRs(VcRS ~ VgRs)^RS _^ ZfgRsVgRS 

dx FgKs dx FgKs Drs 

+ — = 0, (A.9) 

ac«s=— (A.10) 

SgRS=l-ScRS, (A. 11) 

(A. 12) 



PrS — PgRS- 



SgRS^gRsiiRS 

RTgl, 



'■ 0RS 2 
"^wgRS 

where 

FcRS=F,cLi, 
FgRS = Fig + b'is, 

B.C.: 



(A. 13) 



(A. 14) 
(A. 15) 



vfL = vg,s. (A.16) 

In the above, the interfacial friction coefficients are based 
on the following Ergun equation (Tsuo & Gidaspow, 
1990): 



(fig«s4Sc)^(Pc - PgRs). 
. PgRsPc\VgRS — VcRs\ 



EgRsdcSciPc — PgRs) 

\VgRS - VcRs\pcPgRS£cRS 
d,S,{p, - p,Rs) ' 



for EgRs < 0.8, 



for £gRs > 0.8, 
where 

24 



(1 + 0.15/?e°-**'') for Re, 
= 0.44 for Re, > 1000, 

_ \VgRS — ^cRslt^cFgRS 

f^ffRSVgRS^RS 



(A. 18) 
(A. 19) 



The stress modulus of the catalyst is calculated by the 
following equation (Tsuo & Gidaspow, 1990): 
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composition as in Eq. (A.32): 



Gc = 10- 



(A.20) 



The gas-wall friction factor is calculated using the modi- 
fied Hagen-Poisuille equations while the catalyst-wall 
friction factor is calculated using an empirical correlation 
obtained for a commercial scale riser (Kunii & Leven- 
spiel, 1991). 

16 



fgRs = 0.0008 -f 0.0552Ke-° "'' for 10' < Re^ < 10^ 

(A.21) 



{A.31) 
(A.32) 



The temperature dependency of the five kinetic para- 
meters are represented by the following Arrhenius type 
relations: 



= Wexp(-^) .= 1... 



The coke on catalyst is calculated as the sum of the 
residual coke from the regenerator and the coke gener- 
ated during the cracking reactions: 



PgRSyck, 



(A.34) 



where 
Re„ = 



DrS f^gKS PgRS ^gRS 
IJ-gRS 



The weight fractions of the four lumps comprising the 
reaction mixture in the riser can be predicted by the 
following set of component mass balances: 



djnRs 

dx 



Pc£crsOrs<I>c , 



fo, gl, gs, ck, and v 



where 



- (klRS + k2RS + k3Rs)yjoRS, 



= (k.RsyU 



- k4RsyglRS - i^SRSyglRs), 



i^gsRS — i^lRSyjoRS + ^4KS>'sJi<s)> 
Rc,RS = (k^RSyjoRS + l<5RSyg,Rs}, 
R^vRS = 0. 

B.C.: 



(A.26) 
(A.27) 
(A.28) 
(A.29) 



The temperature profiles of the two phases along the 
riser are governed by the following energy balance equa- 
tions: 



t^gRS^pgRS 

+ Pc^cRsQreacRs'}, (A.36) 

= -{hHiRsk.RsyloRS + AH2Rsk2RsyjoRS 
+ AH^Rsl^SRSygoRS + *Rsl<^4^RSyglRS 

+ AHsRsk5RsygiRs)<Pc- (A.37) 



where 

QreacRS 



B.C.: 

Tfj^s = T,Fs and Tf^s = T.^s- (A.38) 

A. 3. Disengaging-stripping section 

The disengaging-stripping section is modeled as a per- 
fectly mixed continuous tank with no reaction taking 
place. In this study, the following exponential type strip- 
ping function is developed to consider the catalyst-to-oil 



The catalyst deactivation function 4)^ in Eq. (A.24) is 
given by Eq. (A.31) where the catalyst deactivation coef- 
ficient is related to the temperature and the feedstock 



r ii.si''cRsi'\.'A 



Then the coke on the catalyst exiting the stripper 
calculated by 



(A.39) 
can be 

(A.40) 
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The mass balances are set up for gases and catalysts in 
the disengaging-stripping section, respectively: 

= •P'jJ!s(l — yfkRs) + Pss — PgRT — PcRsCckST- 

(A.41) 

= F.i«s - Fcsv2 - FcMF- (A.42) 

The weight fractions of the gaseous components exiting 
to the main fractionator are predicted by the following 
equations: 



E 

]=go.gl. and g 



i] 



~ ViRT] F gRT + 



i = 9o,gl and g 



(A.43) 



y^.RT = 1 - ygoRT " VglRT " VgsRT- (A.44) 

The gas oil conversion based on the fresh feed flowing 
into the reactor is given by the following equation: 

V — 1 _ ^gl^^ygoRT . . 



The energy balance is described by the following equa- 
tion: 



{WcRrCpc + 'WgRrCpgRr) — — 

+ FcRsiH'JjiT — ^cRt) + (FgRS — F 
— Fss ^Hstripsrr + QlossRT, 

where 

Ha. = Cp^„(T<fe' - T,,r), 
ti%RT = Cp^giTfi^^ - r,ef), 



HTr't = Cp,{Tj,r - T,.f), 
HZ'rt = CpUTrt - T,.f), 

QlossRT = UnrARTeiTe - T p. 



"hgRT) 
(A.46) 

(A.47) 
(A.48) 
{A.49) 
(A. 50) 
(A. 51) 
(A.52) 
(A.53) 
(A. 54) 



The gas-phase pressure in the disengaging-stripping 
section can be calculated by the following equation of 
state: 



M, 



wgRT 'gRT 



where 

VgRT = Vrt-- 



(A.56) 



The pressure at the bottom of the disengaging-stripping 
section is higher than the pressure at the top by the static 
head exerted by the catalyst: 



A. 4. Reactor cyclones 

The reactor cyclone is modeled as a continuous stirred 
tank (CST). After applying a pseudo-steady state as- 
sumption to the catalyst balance in the reactor cyclones, 
the rate of accumulation of catalyst in the cyclones can be 
expressed by 

IcyiVcCYi^cYii-n KcYi\ =0. 

Vrt L Pt'cri J 

(A.58) 

Then, rearranging the above equation yields the follow- 
ing equation: 

^CYirncYlVcCYlWcRT^CYu V ^ ^^^g^ 
'Id L ^•?CYiVrt J 



where the cyclone inlet velocity is calculated by the 
equation of Rosin, Rammler, and Intelmann (Schweitzer 
et al., 1988): 



^t^gRT^CYli 



dL^^tCYl(pc - WgKT/ygRT)hci 



(A. 60) 



Then the mass flow rates of catalyst particles and gases to 
the main fractionator are calculated by the following 
equations, respectively: 

FcMF = (1 - ricYiXcYi^CYU^^CY,, (A.61) 



,f^Prt - Pmp 
fir(l - ywvRr)- 



(A.62) 
(A.63) 
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A. 5. Overall regenerator 

The equations associated with the total catalyst and 
gas holdups, superficial velocity, and pressure distribu- 
tion are derived on the basis of the overall regenerator 
dynamics. The catalyst and gas holdups in the entire 
regenerator can be calculated by the following equations: 



— t" aD -T r cChl ^ckCLl - 

where 

PgSG = KsaL{XpSG)y/i^ " i 



(A.64) 



(A.65) 



(A.66) 



Since the mass influent to the gas phase increases in 
proportion to the coke burned off the catalyst, the super- 
ficial gas velocity can be correlated to the average gas 
density in the regenerator: 



^ Fgoll - dgRo) + I-'gSG^gRG 
PgRG^RG 



(A.67) 



The average regenerator pressure is assumed to have the 
pressure value at the dense bed exit and can be calculated 
by the following equation of state: 



_ PgRG^Tj}ZgiiG 



= 02,CO,C02,H2 0,N2 



Prc 
where 

M^gRO = X My,ifij, 

and j = D, F, 
fij = C-.jlY^ C^j, i = 02,CO.C02,H2 0,N2 
and j=D,F. 



(A.70) 



(A.71) 



The pressure acting on the bottom of the regenerator can 
be calculated using the static head exerted by the catalyst: 



P«Gb = Prg + -^I^cdHd + jj^cFiz) dz|. (A. 



73) 



A.6. Dense bed 



The dense bed is modeled as a hybrid reactor that 
employs a mixed-tank model for energy and coke bal- 
ances but a tubular reactor model for gas component 
balances. The component balance equations for each 
gaseous phase are described by the following partial 
differential equations, respectively: 

For the emulsion phase: 



dt 



dCi 



iE 



dz 



£gE 



I = 02,CO,C02,H20,N2, 
where 



{CiB — Cie) + RiE, 

(A.74) 



Ro.E 



^1(0.5 + 0.25q)r,E/M^ck 

^gE 

+ (1 + 0.25q)r2E/M„,, + O.Sr^fi] 
■■^LriE/M^,,-r^E'\-r,E. 



Pc£cDq 



['•2£/M„c* + T^.e'] + • 3£;> 



B.C.: 



= 02,CO,C02,H20,N2, 



(A.75) 
(A.76) 

(A.77) 

(A.78) 
(A. 79) 
(A. 80) 

(A.81) 



JfE'PgRG 



= O2, CO, C02,H20,N2, (A.82) 
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For the bubble phase: 






dCiB dCiB CiB^gB 






I = O2, CO,C02,H20,N2 




(A. 83) 












(A.84) 


^COfi = ~ f3B> 




(A.85) 


^COjB — 3B> 




(A. 86) 


^HjOB = 0> 




(A. 87) 


K-N,B = 0, 




(A. 88) 


^ dSgB , d{VgB£0B) 
dt dz ' 




(A.89) 



I.e.: 

Ci^'* = CI^b', i = 02,CO,C02,H20,N2, (A.90) 
B.C.: 

Cf^^=-^r-> i = 02,CO,C02,H20.N2. (A.91) 

In the above equations, the kinetic rates are given by the 
following expressions: 



s = ksRoCcoBC^'BC^, 



(A.92) 

(A93) 
(A.94) 
(A.95) 
(A.96) 



where the temperature dependency of the rate constants 
are assumed to follow the Arrhenius type relations: 

= fciORoexp(^-:^^). =1,3,4 and j = D,F. 

(A.97) 

Once the molar composition of each phase is obtained, 
the mean molar concentration profiles along the dense 
bed can be calculated by taking a weighted average of 
each phase values; 

„ CjEVgE^gE CjBVgBEgB 

'^iD — > 

^gRG 

( = 02,CO,C02,H2 0,N2. (A.98) 



Since the component balance equations presented 
above include the volume fractions and velocities of the 
dense bed phases, they have to be correlated to other 
process variables. First, we explicitly state the require- 
ment that the volume fractions of the catalyst, emulsion 
gas, and bubble gas sum up to unity at any position: 



i + e^B = 1 for 0 < z < Hb. 



(A.99) 



It is evident from the perfect-mixing assumption for the 
catalyst that it has a constant volume fraction. An em- 
pirical correlation between the volume fraction and the 
average superficial gas velocity was developed by a non- 
linear regression of the experimental data of Schnitzlen 
(1987). 

e,D = 0.3418 exp[ - 0.975 Iti^jjc] + 0.1592, (A. 100) 

(A.101) 



1 f"' 

UgRo = -TT- ^bRg{z) dz. 

Jo 



Next, the volume fraction of the bubble gas is predicted 
by the following expression proposed for fast bubbles 
and Geldart A type particles (Davidson & Harrison, 
1963; Kunii & Levenspiel, 1991): 



(A. 102) 



where the minimum fluidization and bubble rising vel- 
ocities are given by 



dcipc - PgRG)g BgmSf 



(A. 103) 
(A. 104) 



dt in Eqs. (A. 103) and (A. 104) denotes the bubble dia- 
meter and is correlated to the axial position in the bed 
(Kunii & Levenspiel, 1991): 



dt = dt,„ - (dt^ 



- d6o)exp( - 0.3z/DnGe), 



(A. 105) 



where the initial bubble diameter dto at the bottom and 
the maximum bubble diameter d^^ at the top of the bed 
are given for a porous plate air distributor as follows 
(Kunii & Levenspiel, 1991): 

dto = — [Cg - v,„-]\ (A.106) 
9 

dtm = 0.59[u<5?j - Vg^r-^D'kL- (A. 107) 

Finally, the interstitial emulsion gas velocity can be easily 
calculated from the mass balance in the bed: 



The holdups of the catalyst and gas in the dense bed 
are given by the following expressions 

WcD = W,KG - WcF - W,cy2> (A. 109) 

= HoQRoii - eccKjjG, (A.l 10) 

from which the height of dense bed can be calculated by 
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(A.lll) 



The uniformity assumption of the catalyst in the bed 
leads to the lumped equation for the coke deposited on 
catalyst as follows: 



where 

FgD = PgRG^RaUgRG. (A.l 13) 

Recalling further that all the substances in the bed are in 
thermal equilibrium, the energy balance also takes the 
form of the following mixed-tank dynamics: 



(w 



+ + QreactD, 

where 

H'gD = CpaiTao — Tref), 
HgD — CpgoiTi) — Tref), 



FcCLziHll - TO + F,n 



(A.l 15) 
(A.l 16) 
(A. 11 7) 
(A.l 18) 
(A. 11 9) 
(A. 120) 



QreacD 



= 02,CO,C02,H20. 



(A.121) 



The molar concentrations of gaseous substances in the 
freeboard are described by the following component 
mass balances: 



dt Egp 8z 



I-tzf) 



PgRG^gF 

i = 02,C0, CO2, H2O, N2, 
where 



(A. 122) 



Ro^F 



-[(0.5 + 0.25q)rir/M^,fc 

Bgp 

+ (1 + 0.25q)r2p/M^,, + 0.5r4F] - O.Sr^p, 

(A. 123) 



Sgp 



Pc^c 



BgpM 

Rn,p = 0, 
I.e.: 

c{?-^> = c%, 

B.C.: 

rC-H-) — r^H'-y 



lriF/My,cu - ] - r3p, 
[r2F/M^,n + r^.p'] + r^p, 
iO.SriF + O.Sripl 



: = 02,CO,C02,H20,N2, 



i = 02,CO,C02,H20,N2. 



(A. 124) 

(A. 125) 

(A. 126) 
(A. 127) 

(A. 128) 

(A. 129) 



In the above equations, the kinetic rates are given by the 
following expressions: 



^IRG^ckF^ 



(A. 130) 

(A. 131) 
(A. 132) 
(A. 13 3) 



The vertical distributions of the catalyst and gas vol- 
ume fractions are estimated as follows. We start again 
with the requirement that the volume fractions of the 
catalyst and gas phases sum up to unity at any position: 



1 + ff ' 
fiF = ripCT, 

r^p = k^RaCcQpCo^pC^^op, 
= /c^jjgCcofCojF- 



= 1 for Ho 



(A. 134) 



The catalyst volume fraction in the freeboard tends to 
decay exponentially with increasing height, and can be 
estimated with the following empirical equation (Li 
& Kwauk, 1980): 
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where 



(A.135) 



-T-f), 



(A. 146) 



(A. 136) 



Then the holdup of the catalysts in the freeboard are 
given by 



4: 



J 

i = cfc,02,C0, COj.HzO, 



(A. 138) n"'^' = Tfo, 



(A. 147) 



(A. 149) 



hno=Ho+Hp. (A. 139) 

The catalyst flow rate in the freeboard is assumed to be 
constant, leading to the following expressions for the 
catalyst flow rate and the superficial catalyst velocity: 



u - 

Pc^RG 



(A. 140) 
(A.141) 



Unlike in the dense bed, the coke on catalyst in the 
freeboard is expressed by a distributed parameter model 
as follows: 



I.e.: 

CSi?> = Cifpo, (A.143) 
B.C.: 

C'if "> = C,,D. (A.144) 

Neglecting the heat transfer resistance between the 
catalyst and gas phases, the temperature distribution in 
the freeboard region can be calculated by the following 
partial differential equation: 



dTp 



(PcScpCpc + PgRGf^gpCpgp)- 



= Qus. 



y- QreacF, 

(A.145) 



J-i.t.HD) 



A. 8. Regenerator cyclones 

The regenerator cyclone is modeled as a continuous 
stirred tank like the reactor cyclone. Similarly to the 
reactor cyclones, the rate of accumulation of catalyst is 
given by 



_ PbCY: 



■Y2d L^''CI'2 ^Cr2 J 



where 

9//gj;G^3c] 



'cr2(-''cy2j 



dlnT^N tCYtiPc — PgRG)hcY7 



(A. 151) 



(A. 152) 
(A.153) 



The flow rate of the catalyst entering the stack through 
the cyclone is given by 



= (1 - riCYlWcCYlhcYl- 



(A. 154) 

The temperature rising across the regenerator cyclone is 
estimated using the following simple correlation 
(Kurihara, 1967; Hovd & Skogestad, 1993): 

TcYi-Tf-' + CTf'&TF- (A. 155) 



A.9. Catalyst transport lines 



The flow rate of catalyst through a catalyst transport 
line is determined by the pressure drop across a slide 
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valve. The mass flow rate of the regenerated catalyst 
through the slide valve can be calculated by 

i'csvi — 

Ksvifu(x,svi yPRGb + Ml - e,cLt)9lcLi sin0cti/lOOO - Pes- 

(A. 156) 

Since the catalyst bulk density is assumed constant in the 
catalyst transport line, the flow rate through the trans- 
port line is just given by 

FcCLi=Fcsvi- (A.157) 

Ignoring the transport lag throughout the catalyst trans- 
port lines, the temperature and the coke on catalyst at the 
outlet of the regenerated catalyst transport line are sim- 
ply given by 

TcCLi = To, (A.158) 

Cc,cLi=C,,u. (A. 159) 

Similar equations can be derived for the spent catalyst 
transport line from the reactor to the regenerator as 



follows: 








y/PRTb + PbCLiglcLl sinecZ.2/1000 - Pro, 


(A. 160) 


FcCL2 — FcSV2! 


(A.161) 


TcCLi — Tut, 


(A.162) 


^ckCL2 — CckRT- 


(A.163) 
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Abstract - 

A dynamic simulator was developed which implements the detailed dynamic model for an FCC process and the model solver- 
presented in Part I of this paper. The simulator incorporates the correlation equations developed in this study for the thermodynamic 
properties and transport parameters contained in our model. First, the simulator was validated by comparing the overall steady-state 
behavior of the system with those in the literature. Then, base case steady-state profiles were obtained for major process variables in 
the reactor riser and regenerator of an industrial scale FCC unit Next, the issue of multiple steady states in FCC operations was 
addressed and confirmed by investigating the sensitivity of our model to initial conditions. Finally, the dynamic responses to step 
changes in three major process inputs were presented and discussed with an emphasis on the interaction between the reactor and 
regenerator dynamics. © 2001 Pubhshed by Elsevier Science Ltd. 

Keywords: Fluid catalytic cracking; Dynamic simulator; Property estimation; Multiple steady states; Sensitivity 



1. Introduction 

Dynamic simulations have been widely used in chem- 
ical process engineering to predict and analyze the dy- 
namic behavior of chemical processes. An accurate and 
efficient dynamic simulator is recognized as a valuable 
tool for various studies on advanced control, online op- 
timization, hazard analysis, and operator training. Since 
a typical FCC unit processes a large amount of feedstock 
into products of additional value, a refinery can achieve 
enormous economic benefits from such process improve- 
ment studies. 

A detailed dynamic model of an FCC process was 
developed in Part I of this paper on the basis of conserva- 
tion principles. Then an efficient model solver was con- 
structed on the basis of a modular approach where the 
model equations were grouped into 12 modules each 
corresponding to a specific part of the process and 
type of equations. The model and numerical algorithms 
developed in the previous part are implemented by a 
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dynamic simulator written in Fortran. Since our dynamic 
simulator for fluid catalytic crackers (called DSim-FCC) 
takes the structure of Fortran 90 module procedures 
comprising versatile subprogram units (Meissner, 1995), 
it can be easily modified when one wants to extend the 
simulator function or to append additional modules. 
Interested readers can obtain the dynamic simulator by 
writing to the authors. 

This paper presents the simulation and its results ob- 
tained using DSim-FCC. First, the estimation technique 
is presented for the thermodynamic properties and 
transport parameters contained in our model. The 
correlation equations for each property estimated on 
the basis of data in the various literatures are presented 
in detail in the Appendix A. Next, to show the vaKdity 
of our sunulator, the steady-state behavior of the 
system is compared with those in the Hterature. Then, 
steady-state profiles of major process variables along 
the axial position of the riser and regenerator are 
shown and compared with other information in the Kter- 
ature for an industrial scale FCC unit. The issue 
of multiple steady states in FCC operations is also 
addressed and confirmed by investigating the sensitivity 
of our model to initial conditions. Finally, the dynamic 
responses to step changes in major process inputs are 
presented and discussed. 
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2. Property estimation 

The dynamic model upon which our simulator DSim- 
FCC is based contains numerous property parameters 
that have to be specified before carrying out simulation 
runs. The parameters mainly consist of the thermo- 
dynamic properties of the substances that are involved 
either in cracking reactions or in coke burning reactions 
and the transport parameters between two phases in 
the reactor riser and the regenerator. Since the 
cracking reaction mixture comprising numerous hydro- 
carbon species is represented by the four lumps (gas 
oil, gasoline, light gases, and coke) in this study, it 
is necessary to estimate the properties of each lump 
using an appropriate characterization method for 
petroleum fractions. It is also necessary to develop 
correlation equations that express the properties as 
functions of system temperature and pressure. This 
section briefly describes the scheme used in this paper 
to develop such correlation. The correlation equations 
incorporated in our simulator are presented in detail in 
the Appendix A. 

True boiling point (TBP) or American Society for 
Testing and Materials D86 (ASTM D86); ASTM Com- 
mittee D-2, 1991 distillations are two standard methods 
to characterize crude oil or petroleum fractions that have 
a wide range of boiling points. Once either distillation 
curve is obtained experimentally, it can be converted into 
the other and can be used to calculate average boiling 
points of various types. The average boiling point thus 
obtained and the specific gravity measured in the lab 
together form the basis for estimating other physical 
properties of the petroleum fraction such as average 
molecular weight, heat capacity, vapor pressure, and so 
on. Hence the physical properties of the gas oil and 
gasoline lumps that are included in our FCC model are 
correlated with the average boiling points and specific 
gravity as well as with the process temperature and 
pressure. It should be noted that the correlation equa- 
tions presented in Appendix A include several empirical 
equations that were developed in this study by nonlinear 
regression of literature data. 

The heat capacities and viscosities of the light gases 
(Ci-C4) and the products of coke burning reactions are 
expressed as functions of temperature only. Since the 
light gases are treated as a single lump in our model, its 
property is assumed independent of composition. For the 
purpose of property estimation, however, it is assumed 
that the light gas lump has typical composition as found 
in field operations and that its property can be calculated 
as a weighted sum of component properties. 

Our dynamic model postulates two-phase tubular re- 
actor models for the reactor riser and as well as the dense 
bed and freeboard of the regenerator. The parameters 
associated with the heat and mass transfer between these 
two phases are thermal conductivities, heat transfer coef- 



ficients, diffusivities, and mass transfer coefficients. As 
shown in Appendix A, these parameters are estimated 
using correlation equations either found in the literature 
(Baird & Rice, 1975; Kunii & Levenspiel, 1991) or ob- 
tained by nonlinear regression of literature data (Tech- 
nical Data Committee, 1988). 



3. Simulations and discussions 

The simulator DSim-FCC was used to investigate the 
behavior of a side-by-side type FCC unit shown in Fig. 1. 
Recall that the dynamic model presented in Part I of this 
paper was found to have 13 degrees of freedom that has 
to be specified before carrying out an open-loop simula- 
tion. Table 1 lists 13 variables chosen in this study as 
specifications of the degrees of freedom as well as eight 
parameters related to feedstock characteristics and six 
parameters related to ambient air characteristics. From 
the viewpoint of process control, these 27 variables rep- 
resent either external disturbances to the process or pro- 
spective manipulated variables in a certain control 
configuration. From the viewpoint of dynamic simula- 
tion, they represent process inputs whose values can be 
altered during the execution of the simulator in order to 
obtain dynamic responses to such changes. Table 1 also 
lists the assumed values of the 27 inputs at the base case 
operating conditions. Table 2 lists the dimensions 
of an 30,000 BPD FCC unit assumed in this study. 
Table 3 shows the properties of the catalyst and coke, 
and Table 4 lists the kinetic parameters of cracking 
reactions and coke burning reactions. The kinetic para- 
meters of the four-lump cracking reactions in Table 4 
were reported by Ancheyta-Juarez, Lopez-Isunza, 




Fig. 1. Schematic diagram of a side-by-side type FCC unit. 
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Table 1 

Process inputs and base case operating conditions 



Type Variable Description Base case operating 

conditions 



Degrees of freedom 



Feedstock characteristics 



ir characteristics 



Entrance flow rate of liquid feedstock 49.3 kg/s 

Entrance temperature of liquid feedstock 535 K 

Flow rate of the air entering the regenerator 34.0 kg/s 

Temperature of the air entering the regenerator 432 K 

Flow rate of stripping steam 1.1 kg/s 

Temperature of stripping steam 430 K 

Flow rate of dispersion steam 1.6 kg/s 

Temperature of dispersion steam 430 K 

Pressure at main-fractionator 185.5 kPa 

Downstream pressure of the stack gas valve 1 10.0 kPa 
Stem position of the slide valve on the 

regenerated catalyst transport line 0.57 
Stem position of the slide valve on the spent 

catalyst transport bne 0.46 

Stem position of the regenerator stack gas valve 0.70 

TBP distillation temperatures at distilled vol% 
equal to 10, 30, 50, 70, and 90, 
respectively. 

Specific gravity of liquid feedstock 
Aromatics to naphthenes weight ratio in feedstock 
Weight fraction of Conradson carbon residue in feedstock 0.5 wt% 



554.3,605.4, 
647.0, 688.2, 

and 744.8 K, respectively. 
0.894 (API: 26.8) 
2.1 



AV,f(:o\f<&\ Composition of ambient ai 

T, Ambient air temperature 



0.2097, 0, 0.0003, 0 
0.79, respectively. 
300 K 



Riser 30.0 1.1 

Disengaging section 8.5 5.1 

Stripping section 6.2 3.2 

Regenerator 18.0 8.2 

Spent catalyst transport line 10.5 0.8 

Regenerated catalyst transport line 12.5 0.8 

Total no. of cyclones equipped Dimensions of cyclon 

Reactor cyclones 2 0.8 (H)/0.4 (W) 

Regenerator cyclones 2 1.0 (H)/0.5 (W) 



Table 3 

Characteristics of catalyst and coke used in simulations 
Catalyst 

Average particle diameter 7 x 10 ~ ' m 

Minimum particle diameter 2 x 10 " ' m 

Sphericity 
Density 
Specific heat 

Coke 
Density 
Specific heat 

Atomic ratio of hydrogen to carbon 



0.72 

1410 kg/m^ 
1.15kJ/(kgK) 



IS catalyst 
IS catalyst 



Aguila-Rodriguez, and Moreno-Mayorga (1997). The 
heats of cracking reactions were roughly estimated by 
trial-and-error to match the calculated riser outlet tem- 
perature with a typical operation temperature reported 
by CaUahan and Ushiba (1988). The kinetic parameters 
for oxidation of coke, homogeneous and catalytic oxida- 
tions of carbon monoxide were adopted from the papers 
of Morley and de Lasa (1987), Howard, Williams, and 
Fine (1973), and Ali and Rohani (1997), respectively. The 
intrinsic CO2/CO molar ratio was estimated using the 
empirical equations proposed by Errazu, de Lasa, and 
Sarti (1979). 
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Table 4 

Kinetic parameters of cracking and coke burning 



Four-Lump cracking n 



Gas oil to 
Gas oil to 
Gasoline 1 



C1-C4 gases 

o C1-C4 gases 
o coke 



Frequency factor (s ^) 

1457.50 

127.59 



Catalyst deactivation 
Coke burning reaction 



Activation energy (kJ/kg m 

57359 

52754 

31820 

65 733 

66570 

£1. = 49000kJ/kgmol 



and oxidation 



Coke burning r« 
CO oxidation (non-catalytic) 
CO oxidation (catalytic) 
Intrinsic COj/CO molar 



Heat of formation 



Frequency factor 
1.4x 10*m7(kgmol s) 
3.5x 10'm7(kgmols) 
247.75 m* V(kgmol° ' kg s) 



Activatii 
125 000 
165000 
70480 



n energy (kJ/(kg mol) 



a = 0.000953 exp(5585/r) 

AHj.,» = - 4800.22 + 16.1 T, kJ/kg mol 

AHyo. = - 10 364.88 + 34.60 T + 0.00055 kJ/kgmol 

AH/co = - 118975.04 + 27.6ir + 0.00251^^ kJ/kgmol 

AH/CO. = - 406909.11 + 43.26r + 0.0057ST^ kJ/kgmol 

AH^.o = - 252111.38 + 34.39r + 0.000315T^ kJ/kgmol 



The simulation results presented in this section cover 
both the steady state and dynamic responses of the FCC 
unit described previously. First, steady-state responses 
were obtained by executing the simulator for a suffi- 
ciently long time until aU the variables lined out. These 
steady-state simulations were carried out many times 
each with altered values of the two major operating 
variables (the catalyst circulation rate (f c) and the air 
flow rate into the regenerator Then, the simulation 
results were compared with those in the literature. In 
addition, the predicted profiles of several process vari- 
ables are plotted along the height of reactor riser or 
regenerator for the base case operating conditions. Next, 
the multiple steady states in FCC operations were con- 
firmed by investigating the sensitivity of our model to 
initial conditions. Finally, the dynamic responses of the 
process to step changes in several process input variables 
were predicted and analyzed. 

3. 1. Comparison of steady-state behavior with the previous 
models 

In this section we show the validity of our simulator 
by comparing our simulation results with those 
from the previous FCC models. Although direct and 
qualitative comparison is not feasible due to different 
kinetics and unavailability of simulation data, we 
investigate the overall steady-state behavior of the FCC 
unit as functions of the catalyst circulation and the air 
flow rate. 

Fig. 2 shows the steady-state behavior of major 
state variables as functions of the air-flow rate when the 



catalyst circulation rate is fixed at a specified value of 250 
or 300 kg/s. Fig. 2a shows that with decreasing air flow 
rate the oxygen concentration gradually decreases until 
almost aU the oxygen is exhausted. In this transition the 
system moves from a so-called full-combustion mode to 
a partial-combustion mode (Arbel, Huang, Rinard, Shin- 
nar, & Sapre, 1995a; Arbel, Rinard, Shinnar, & Sapre, 
1995b). The point of transition in combustion mode can 
be easily located at the maximum points of the curves for 
the regenerator and reactor temperatures in Fig. 2b. 
Alternatively, this point corresponds to either the max- 
imum points of the gas-oil conversion or the minimum 
points of the coke on catalyst in Fig. 2c. The results 
shown in Fig. 2 exhibit good matches with the field 
observations and simulation results of Arbel et al. (1995a) 
not only in terms of qualitative trend but also in terms of 
the slopes of the curves. 

Fig. 3 shows the steady-state behavior of major state 
variables as functions of the catalyst circulation rate 
when the air flow is fixed at a specified value of 35 or 
40 kg/s. As the catalyst circulation rate increases, the CO 
combustion in the regenerator transits from a full com- 
bustion mode to a partial combustion model with ac- 
companying changes in the state variables. Again, these 
results are in good agreement with Arbel et al. (1995a). 
The steady-state behavior of the regenerator temperature 
shown in Fig. 3b is also analogous to that of Kumar, 
Chadha, Gupta, and Sharma (1995). Ali and Rohani 
(1997) presented similar results for the regenerator tem- 
perature and the coke on catalyst, but opposite trend in 
the reactor temperature to ours as well as Arbel et al. 
{1995a). 




Fig. 3. Steady-state responses to changes in catalyst circulation rate. 
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3.2. Steady-state profiles along the riser 
and the regenerator 

Fig. 4 depicts the predicted steady-state profiles of 
several process variables along the reactor riser under the 
base case operating conditions. The conversion of gas oil 
reaches 73 wt% at the riser outlet and 90% of the conver- 
sion is attained within the first 20 m of the riser (Fig. 4a). 
Most heat exchange between the catalyst phase and the 
gas phase takes place within the first 10 m, thus two- 
phase temperatures approaching each other within 3.0 K 
at the riser outlet (Fig. 4b). After sufficient heat is trans- 
ferred from the catalyst to the gas phase, the gas temper- 
ature as well as the catalyst temperature drops with 
increasing riser height because of endothermic cracking 
reactions. Both the catalyst and gas velocities sharply rise 
to about 7.5 m/s at the riser bottom due to the volume 
expansion accompanying feed vaporization (Fig. 4c). 
There is significant molar expansion of gaseous sub- 
stance as large molecules crack to smaller molecules, so 
the interstitial velocities of both phases steadily increase 
up to 17 m/s at the outlet. The slip velocity between two 
phases is maintained within 0.25 m/s. It can be inferred 
from the predicted velocity profiles that assuming con- 
stant velocity in the riser would result in significant error 
in describing cracking reactions. The pressure shows an 
almost linear decrease along the riser with the total drop 
of about 1 6 kPa. This value shows good agreement with 



a measured value (16.7 kPa) and a predicted value 
(13.7 kPa) reported by Theologos, Nikou, Lygeros, and 
Markatos (1997) who developed a very rigorous three- 
dimensional riser model. 

Fig. 5 shows the predicted profiles of several process 
variables along the regenerator height at the same steady 
state as in Fig. 4. Since the predicted catalyst bulk density 
of 448 kg/m^ in the dense bed is much larger than 
3 1 kg/m^ in the freeboard, most coke burning reactions 
take place in the dense bed (Fig. 5a). The molar stack-gas 
composition was predicted to be 0.2% O2, 0.6% CO, 
14.8% CO2, and 9.2% H2O, which also agrees well 
with the observations obtained for typical promoted 
coke burning reactions (Faltsi-Saravelou, Vasalos, 
& Dimogiorgas, 1991; Zheng, 1994). Since catalyst par- 
ticles are assumed perfectly mixed in the dense bed, they 
have uniform coke concentration and temperature pro- 
files in the dense bed (Figs. 5b and c). As a result of the 
after-burning oxidation of CO, CO concentration is 
lower and temperature is higher in the freeboard than in 
the dense bed. 

3.3. Multiple steady-states and sensitivity to 
initial conditions 

There has been plethora of literature on the subject of 
multiple steady states and their stability of FCC units in 
the past three decades. Since Iscol (1970) proposed the 
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existence of multiple steady states in the normal ranges of 
commercial operation, many papers followed arguing 
pro and con on the possibility of multiplicity and its 
conditions (Lee & Kugelman, 1973; Edwards «fc Kim, 
1988; Elshishini & Elnashaie, 1990, Arandes & de Lasa, 
1992). The dispute and confusion on the issue seems to be 
clarified by Arbel et al. (1995a,b), who performed a sys- 
tematic analysis based on a relatively detailed model of 
modern FCC units. A part of the authors' conclusions 
can be summarized as follows: (1) there are at least three 
steady states for fixed inputs to present design units. In 
case there is only one, it is a trivial quenched state; (2) the 
upper (ignited) and the lower (quenched) state are always 
stable whereas the intermediate one is always unstable. 

Most analyses on multiple steady states for a process 
involving exothermic reactions are carried out using 
a steady-state model to find the operating points at which 
the rate of heat generation equals that of heat removal. 
Dynamic models can be also used to confirm the multiple 
steady states by investigating the trajectories of the sys- 
tem subject to different initial conditions because each 
stable steady state has a unique domain of attraction on 
the phase plane. Elnashaie and Elshishini (1993) studied 
the sensitivity of their dynamic model to initial condi- 
tions. The process trajectories starting from several initial 
conditions were found to drift either to an upper steady 
state or to a lower steady state. They found that the 
dynamic behavior of their FCC model is very sensitive to 
the initial catalyst activity in the reactor. 



Dynamic simulations were performed starting from 
perturbed initial conditions to investigate multiple steady 
states for our FCC model. First, the initial conditions 
were set as the values at the base case steady state 
obtained previously except the dense bed temperature 
which was perturbed to 1 100, 900, and 820 K, respective- 
ly. Fig. 6 compares the transient trajectories of reactor 
temperature, dense-bed temperature, and coke on regen- 
erated catalyst evolving from each initial condition. 
When starting from 1100 or 900 K, the process shows 
different transient responses (curves 1 and 2) but ulti- 
mately converges to the original base case steady state. 
When starting from 820 K, however, both the reactor 
and dense-bed temperatures keep decreasing until it 
reaches the point where the regenerated catalyst is not 
capable of vaporizing gas oil feed any more. Since then 
the riser wiU be filled with liquid gas oil feed, normal 
operation of the unit would not be possible any longer. 
This state of inoperability of an FCC unit is called 
a quenched state in the literature (Elshishini & Elnashaie, 
1990; Arbel et al., 1995b) whereas the normal steady state 
obtained above is called an ignited state. It is clear that 
the unstable steady state in the middle cannot be identi- 
fied in open-loop dynamic simulations. 

It may be interesting to check the sensitivity of 
our FCC model to initial conditions of several key vari- 
ables. It was found in our simulations that the process is 
sensitive to dense bed temperature in that an initial 
temperature below 827 K would cause the process to 
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Fig. 6. Transient responses starting from different initial dense bed temperatures: (1) Tt°' = 1100 K, (2) rt°> = 900 K. (3) Tj?' = 820 K. 



drift toward the quenched state. On the other hand, the 
process always drifts toward the upper steady state what- 
ever temperature the reactor initially may assume be- 
tween 270 and 920 K. The different sensitivity behavior 
between the two cases can be attributed to much larger 
thermal holdup of the regenerator compared to the reac- 
tor. The process also drifts to the upper steady states 
when the initial coke on spent and regenerated catalyst is 
less than 0.055 and 0.054, respectively. 

It is the regenerator with exothermic reactions and 
back-mixed thermal dynamics that is responsible for the 
multiple steady states in FCC operation. Since the riser is 
modeled as a tubular reactor with endothermic reactions 
and the disengaging-stripping section does not involve 
any reaction, the reactor does not have dynamic charac- 
teristics that may lead to multiple steady states (Maya- 
Yescas & Lopez-Isunza, 1997). To investigate the 
steady-state multiplicity inherent in the regenerator dy- 
namics, dynamic simulations were performed with the 
scope confined to the regenerator. Fig. 7 compares the 
transient responses evolving from two different initial 
temperatures (1 100 and 500 K) of the dense bed when the 
conditions of spent catalyst entering the regenerator are 
fixed. It is clear that the high-initial temperature leads 
to the upper steady state whereas the low initial 
temperature leads to the lower steady state. The upper 
and lower steady states represent two stable modes of 



Ltppsr steatly-^stafe 



Lower steady-stafe 



r 



ijpper steady-state : 



1 regenerator operation with catalyst 



I.-S. Han, C.-B. Chung I Chemical Engineering Science 56 (2001) 1973-1990 



Table 5 

Step change scheme for process inputs in dynamic simulation 



Case 1: Changes in gas oil feed rate. 
Case 2: Changes in air flow rate. 

Case 3: Changes in stem positions of the slide valves at both catalyst transport lines. 



operation as far as the regenerator is concerned. When 
connected to a reactor with catalyst circulated in-be- 
tween, however, the lower steady state would result in the 
quenching of the whole unit because the regenerated 
catalyst at such low temperature would fail to vaporize 
gas oil feed. 

3.4. Dynamic responses 

Dynamic simulation of the FCC process was per- 
formed according to the simulation scheme shown in 
Table 5. Gas oil feed rate, air flow rate, and stem posi- 
tions of the slide valves at both catalyst transport lines 
were chosen as simulation variables to whose changes the 
dynamic responses are demonstrated. Each simulation 
run started from the steady state corresponding to the 
base case operating conditions and the subsequent tran- 
sient response was obtained as each simulation variable 
went through a series of step changes shown in Table 5. 



3.4.1. Changes in gas-oil feed rate 

Fig. 8 shows the predicted dynamic responses to the 
changes in gas oil feed rate. After the gas oil feed rate is 
increased by 5% at time equal to lOmin, the reactor 
temperature drops because of an increase in heat con- 
sumption to vaporize the gas-oil feed (Fig. 8a). The 
lowered riser temperature results in a considerable de- 
crease in gas oil conversion by about 6.1 wt% (Fig. 8e) 
while the increased feed rate leads to a sharp rise in coke 
on spent catalyst and subsequently on regenerated cata- 
lyst (Fig. 8f ). Because the increased feed rate also raises 
the density of hydrocarbon gases in the reactor, there is 
a jump of gas pressure as well as of the pressure at the 
reactor bottom (Fig. 8c). This jump will also trigger 
a jump of spent catalyst flow rate but an opposite effect 
upon regenerated catalyst flow rate (Fig. 8d). Therefore, 
the catalyst holdup in the reactor steadily decreases until 
a new steady state is reached where the two catalyst 
circulation rates are balanced each other (Fig. 8b). Since 
decreasing catalyst holdup means decreasing static head 



acting on the reactor bottom, the pressure at the bottom 
will descend from the aforementioned initial jump until it 
settles down to a new steady-state value which is lower 
than the original steady state (Fig. 8c). This inverse re- 
sponse behavior of the reactor bottom pressure also 
explains the similar behavior of the spent catalyst flow 
rate which is directly influenced by the reactor bottom 
pressure (Fig. 8d). Regenerator temperatures start to de- 
crease in both the dense bed and freeboard regions since 
the increased spent catalyst inflow from the reactor 
cooler than before implies an enhanced heat removal rate 
in energy balance (Fig. 8a). Since the fixed air flow does 
not supply enough oxygen to burn ofl^ the increased coke 
in the dense bed nor to sustain after -burning reactions in 
the freeboard, CO concentration in the stack gas steadily 
rises (Fig. 8e). The oxygen deficiency in the regenerator 
also leads to no more rise in the freeboard temperature 
above the dense bed temperature at the new steady state 
(Fig. 8a). 

When the gas oil feed rate is decreased to its previous 
value at time equal to 150 min, the system completely 
recovers its original steady state. The transient responses 
over the second time interval in Fig. 8 show in what an 
orderly and coordinated fashion the process variables are 
returning to the old steady state. 

A 5% step decrease in the gas oil feed rate at 300 min 
causes the system to move in the opposite direction to the 
case of the 5% step increase but the detailed responses 
show quite different aspects due to process nonlinearity. 
The gas oil conversion increases by about 3.9 wt% due to 
increased catalyst-to-oil ratio in the riser as well as better 
regeneration of catalyst in the regenerator (Figs. 8d-f). 
The reactor bottom pressure and the spent catalyst flow 
rate again show inverse responses which, on this occa- 
sion, are characterized by initial drop followed by event- 
ual rise (Figs. 8c and d). Accordingly, catalyst holdup 
shows a steady increase in the reactor whereas a 
steady decrease in the regenerator (Fig. 8b). Contrary 
to the oxygen deficiency encountered before, oxygen 
surplus prevails throughout the regenerator because of 



1982 



I.-S. Han, C.-B. Chung / Chemical Engineering Science 56 (2001) 1973-1990 




Reactor bcabm 





(a) 


» Spent 






Fig. 8. Dynamic responses to step changes in gas oil feed rate (case 1). 



the decreased spent catalyst inflow from the reactor. This 
causes a rise in the dense bed temperature, a much higher 
rise in the freeboard temperature due to after-burn of 
CO, and a consequent drop in the CO concentration of 
stack gas (Figs. 8a and 8e). 

It was found in the simulations described above that 
the step changes in gas oil feed rate immediately alter the 
reactor pressure which subsequently triggers a series of 
changes upon the process dynamics. Since the stem posi- 
tions of the slide valves on catalyst transport lines are 
held fixed in our open-loop simulation, catalyst circula- 
tion rates are solely determined by pressure differences 
across the valves. Accordingly, the pressure variation 
leads to serious unbalance between the flow rates of spent 
and regenerated catalyst and consequently to significant 
changes in the catalyst holdups in two vessels. In our 
simulations, the process was found to keep on stable 
operation in the face of step changes in feed rate of magni- 



tude of about ± 7%. However, a large change exceeding 
this range caused the process to drift to a quenched state 
where normal operation of the unit is impossible. 

3.4.2. Changes in air flow rate 

Fig. 9 shows the predicted dynamic responses to the 
changes in the air flow rate to the regenerator. A 5% step 
increase in air flow causes a sharp increase in the regener- 
ator pressure as well as in the regenerated catalyst flow 
rate to the riser bottom (Figs. 9c and d). On the other 
hand, the spent catalyst flow rate shows an inverse re- 
sponse behavior of initial decrease followed by ultimate 
increase (Fig. 9d). This behavior can be explained by 
a steady increase in the reactor bottom pressure (Fig. 9c) 
due to increasing static pressure exerted by the rising 
catalyst holdup in the reactor (Fig. 9b). The increased air 
flow rate also accelerates coke burning and thus raises 
temperature in every part of the unit (Fig. 9a) and reduces 
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Fig. 9. Dynamic responses to step changes in air flow rate (case 2). 



the CO concentration in the stack gas (Fig. 9e). It also 
causes the coke concentrations on both spent and regen- 
erated catalyst to steadily drop to their respective steady 
state values (Fig. 9f). The conversion of gas oil is pro- 
moted by about 3.2 wt% by the elevated riser temper- 
ature and cleaner regenerated catalyst (Fig. 9e). When the 
air flow rate is decreased to its previous value, however, 
all the process variables are completely brought back to 
their original steady state. 

A 5% step decrease in the air flow rate at time equal to 
300 min brings about the state of oxygen deficiency in the 
regenerator. A large amount of coke cannot be burnt off, 
resulting in gradual accumulation on catalyst surfaces in 
both vessels over a long time (Fig. 9f, where the responses 
are partially shown). Such a long response time is charac- 
teristics of the composition dynamics in the regenerator 
that has large mass and thermal holdups and strong 
interactions with the reactor. CO concentration also 
shows a considerable increase because of poor after- 



burning reactions in both the dense bed and freeboard 
(Fig. 9e). This oxygen deficiency lowers the temperatures 
of both vessels by about 23-39 K (Fig, 9a) as well as the 
conversion of gas oil by 6.3 wt% (Fig. 9e). It was also 
found in our simulations that the process reaches a new 
stable state in the face of step changes in air flow rate of 
magnitude up to about + 8%. 



3.4.3. Changes in the stem positions of slide valves 

Fig. 10 shows responses due to the changes in both the 
spent and regenerated catalyst flow rates as effected by 
slide valve manipulation. The stem positions of the slide 
valves on both the transport lines are adjusted upward by 
5% at time equal to 10 min which corresponds to 8.5% 
increase in the catalyst circulation rate. Since both the 
spent and regenerated catalyst flow rates are raised by 
the approximately same amount, the catalyst holdup in 
each vessel is maintained constant (Fig. 10b). There is 
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also little variation in the pressures of the reactor and 
regenerator (Fig. 10c). However, the regenerator temper- 
atures steadily drop both in the dense bed and freeboard 
because the increased flow of spent catalyst implies en- 
hanced heat removal in energy balance for the regener- 
ator. On the other hand, the reactor temperature initially 
rises due to increased catalyst flow from the regenerator 
but drops in a while as the dense bed temperature is 
lowered (Fig. 10a). The CO concentration slightly rises 
because after-burning reactions are suppressed by 
lowered freeboard temperature (Fig. lOe). Consequently, 
the cokes on both the regenerated and spent catalyst 
increase to new steady-state values (Fig. lOf). The 
gas oil conversion slightly increases by about 
1.2 wt% due to an increased catalyst-to-oil ratio and 
elevated reactor temperature (Fig. lOe). If the stem 
positions of the slide valves are shifted back to their 
original positions, the process variables converge again 
to their previous values. 



On the other hand, a 5% decrease in the stem 
positions at 300 min reduces the catalyst circulation rate 
by 9.6%. This reduction in the catalyst circulation rate 
decreases the gas oil conversion by 1.1% (Fig. lOe). The 
regenerator temperatures increase because the reduced 
spent catalyst flow provides not only a lowered heat 
removal rate but also a new oxygen-rich environment 
favorable to coke-burning reactions (Figs. 10a and e). 
Accordingly, there is a slight decrease in both the coke on 
regenerated catalyst and the CO concentration in the 
stack gas (Figs. lOe and f). 



4. Conclusions 

A dynamic simulator for fluid catalytic crackers 
called DSim-FCC was developed which implements the 
detailed dynamic model for an FCC process and the 
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model solver presented in Part I of this paper. The 
correlation equations were prepared for the thermodyn- 
amic properties and transport parameters contained in 
our model either by literature survey or by nonhnear 
regression of literature data. 

The steady-state behavior was investigated as func- 
tions of the catalyst circulation rate and the air flow rate 
and the results showed good matches with those in the 
literature. For an industrial scale FCC unit of 30,000 
BPD capacity, the base case steady-state profiles 
were obtained for velocity, pressure, temperature 
and composition in the reactor riser and the regenerator. 
Some steady-state predictions showed satisfactory 
agreement with the literature data. Then the issue of 
multiple steady states in FCC operations was addressed 
and confirmed by investigating the sensitivity of 
our model to initial conditions. Depending on initial 
conditions, the process was found to drift either toward 
an upper (ignited) steady state or toward a lower 
(quenched) steady state where normal operation of the 
unit is impossible. It was also found that the regenerator 
with exothermic reactions and back-mixed thermal dy- 
namics is responsible for multiple steady states in FCC 
operation. 

The dynamic step responses were obtained in an 
open-loop mode for gas oil feed rate, air flow rate, and 
stem positions of the slide valves on catalyst transport 
lines. The presented responses show typical inverse re- 
sponses and nonlinear behavior which stem from com- 
plicated interaction between the reactor and regenerator 
dynamics. 
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Notation 






A, S, C 


coefficients of the Antoine equation 


YV 


API 


API gravity 




C 


mole concentration, kg mol/m^ 


y 




coke on catalyst, kg-coke/kg-catalyst 


r,. 


CckSTO 


minimum coke content attainable by strip- 






ping, kg-coke/kg-catalyst 






heat capacity, kJ/(kg K) 


Greek letters 




mean heat capacity, kJ/(kg K) 




d 


average diameter, m 




D 


diameter, m 




Df 


effective difJusion coefficient, m^/s 




E 


activation energy, kJ/kg mol 






exponent in the stripping function, s/kg 


p. 


/ 


mole fraction 




F 


mass flow rate, kg/s 




9 


acceleration due to gravity, 9.8 m/s^ 


AH„ 


hp 


interface heat transfer coefficient between 






catalyst and gas phases, kJ/(m^ s K) 




k 


heat conductivity, kJ/(s m K) 




K 


valve flow rating factor, kg/(s kPa°'*), 


>lc\i> Icy 2 



flow rating factors of the cyclone dip-legs, 
kg/(sm^') 

Watson characterization factor, K^^' 
interchange coefficients between bubble 
and emulsion phase, 1/s 
molecular weight 

numbers of cyclones in the reactor and 
regenerator, respectively 
pressure, kPa 

pseudo-critical pressure, kPa 

pseudo-reduced pressure, kPa 

aromatics to naphthenes ratio in a liquid 

feedstock 

specific gravity 

Temperature, K 

ASTM D86 distillation temperatures at 
distilled vol% equal to 10,30,50,70, and 90, 
respectively, K 

ASTM D86 distillation temperature, K 

initial dense bed temperature, K 

ambient temperature, K 

molal average boiling temperature, K 

mean average boiling temperature, K 

pseudo-critical temperature, K 

pseudo-reduced temperature 

reference temperature, 298.15 K 

TBP distillation temperature, K 

normal boiling point, K 

volume average boiling temperature, K 

superficial velocity, m/s 

overall heat transfer coefficient between 

a process unit and the surroundings, 

kJ/(m^ s K) 

interstitial velocity, m/s 

superficial gas velocity at a minimum 

fluidizing condition, m/s 

holdup, kg 

valve stem position, [0-1] 
weight fraction 

weight fraction of Conradson carbon resi- 
due in feedstock 



catalyst deactivation coefficient 

(valve head differential at maximum 

flow)/(valve head differential at zero flow), 

[0-1] 

catalyst decay constant, 1/s 
heat of formation, kJ/kg mol 
heat of stripping, kJ/kg 
heat of vaporization, kJ/kg 
volume fraction, [0-1] 
void fraction at a minimum fluidizing con- 
dition 

collection efficiencies of cyclones, [0-1] 
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Subscripts 



CL2 
CO 

C02 

CYl 
CY2 
ds 
D 



O2 
RG 
RT 



ST 
SVl 
SV2 



viscosity, kg/(m s) 
pseudo-critical viscosity, kg/(m s) 
pseudo-reduced viscosity 
density, kg/m^ 

intrinsic CO 2 /CO molar ratio in coke 



bubble or catalyst bulk 
bubble phase 
catalyst 
coke 

regenerated catalyst transport line (to 
reactor) 

spent catalyst transport line (to regener- 
ator) 

carbon monoxide 
carbon dioxide 
reactor cyclone 
regenerator cyclone 
dispersion steam 
dense bed 
emulsion phase 
freeboard 
gas 

gasoline 
gas oil 

C1-C4 (lighter) gases 
water 

liquid feedstock (gas oil) 
mixture 

main-fraction ator 

nitrogen 

oxygen 

regenerator 

reactor 

stack gas valve 

stripping steam 

stripping section 

regenerated catalyst slide valve 

spent catalyst slide valve 



(Riazi, 1988): 
TxBp = o(Tastm)'' 

^ASTTM — ^ '''(^TBp)^'*! 



(A.1) 
(A.2) 



where the coefficients a and b vary with the percent of 
distilled volume as shown in Table 6. Then the ASTM 
D86 distillation curve can be used to estimate the vol- 
ume, molal, and mean average boiling points, respective- 
ly (Technical Data Committee, 1988) 

TvABP = 0.2(Tio + T30 + Tso + T-,0 + ^90), (A.3) 

Tmabp = TvABP - 0.5556 exp[ - 0.5638 

- 0.0080(1. 8 TvABP - 491.67)°-^^" 

+ 3.0473(5/)°""], (A.4) 

TweABP = TvABP " 0.5556 exp[ - 0.9440 

- 0.0087(1. STvABP - 491.67)° "" 

-I- 2.9972(50" ""], (A.5) 

where 

(S/) = 0.0225(T9o - Tio). (A.6) 

The Watson characterization factor is calculated using 
specific gravity and mean average boiling point as fol- 
lows: 



(API) + 131.5' 

_ (1.8TMeABp)^^^ 



(A. 8) 



The molecular weights of gas oil and gasoline lump 
can be estimated by Eq. (A.9) (Riazi, 1988), and the 
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Appendix A. Correlation of physical properties and trans- 
port parameters 

ASTM D86 and TBP distillation curves can be con- 
verted into each other by the following expressions 



Table 6 

Coefficients in Eqs. (A.l) and (A.2) 



Vol% distilled 



0 0.9167 1.0019 

10 0.5277 1.0900 

30 0.7429 1.0425 

50 0.8920 1.0176 

70 0.8705 1.0226 

90 0.9490 1.0110 

95 0.8008 1.0355 
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average molecular weight of n component gas mixture is 
calculated by Eq. (A. 10): 



= 42.965[exp(2.097 x lO-^TMeABP 
- 7.787S, + 2.085 x lO^TMeABpS,)] 

„ t1-26007c'4. 98308 
^ MeABP '^g > 



(A.9) 
(A. 10) 



The liquid- and gas-phase heat capacities of gas oil and 
gasoline lump are estimated using a modified Lee-Kesler 
correlation (Lee and Kesler, 1988) which neglects the 
pressure effect since typical FCC processes are operated 
under relatively low pressures ranging between 200 and 
300 kPa: 



= «! -^«2T-t-a3T^ 
= )5j + PiT + PiT\ 
re 

= - 4.90383 + (0.099319 + 0.104281SJK;. 
(4.814066 -0.194833Ky) 

= (7.53624 X 10-'*)(1.0 + 0.82463Xj.) 
0.27634 \ 



1.12172 - 



J' 



as = ( - 1.356523 x 10-'')(1.0 + 0.82463Ky; 

/ 0.70958\ 
X (^2.9027 ^j, 

Pi = - 1.492343 + 0.124432Ky + 



= ( - 7.53624 X 10-'*)[2.9247 - (1.5524 

f 5.0694 1 
-0.05543K^)K^ + ^4^(6.0283 — ) 



, = (1.356523 X 10"*)(1.6946 -I- 0.0884)54), 



(A. 17) 
(A. 18) 



- = [{^-'»)('»-f7) 

x(S^ - 0.885)(5j - 0.70)(10'')]2 

for 10.0 <Kf < 12.8, 

/?4 = 0 for all other cases. (A. 19) 

The heat capacity of the light gas lump is calculated 
by Eq. (A.20) assuming the composition (H2 = 0.2%, 
CI = 5.7%, C2 = 7.8%, C3 = 29.7%, and C4 = 56,6%) 
found in typical field operations (Mcketta and Cunning- 
ham, 1976): 

Cpg, = 0.2457 -I- 5.3 x lO'T - 2.1527 x lO'^T^ (A.20) 

The heat capacities of carbon monoxide, carbon dioxide, 
nitrogen, oxygen, and steam are calculated by the follow- 





ing: 




(A. 11) 


Cpo, = 1.081 -1- 0.000034T, 


(A.21) 


(A. 12) 


Cpco = 0.986 -f 0,0001 8 r, 


(A.22) 




Cpco, = 0.983 + 0.00026r, 


(A.23) 




CpH.o = 1.91 -f 0.000035 T, 


(A.24) 




Cp^v. = 0.971 + 0.0001 5T. 


(A.25) 


(A. 13) 


Then the mean heat capacity of 


n component mixture can 



be calculated by 

i y<Cp,(r)dT 



Cp„ = 



(A.26) 



The gas phase viscosities of hydrocarbon lumps are 
estimated by Eq. (A.27) which factorizes the viscosity into 
pseudo-reduced viscosity and pseudo-critical viscosity 
(Nelson, 1958) 



= 3.515 X 10- 



(A.27) 



The correlation equation (A.28) for pseudo-reduced vis- 
cosity was obtained by nonlinear regression of literature 
data (Nelson, 1958) over the range of Tp, and Pp, given in 
Eq. (A. 29) 

Hp, = 0.435 exp[(1.33 16 - r°;«*'^')i^p,]Tp, 

+ 0.0155, (A.28) 

0.75 < Tp, = ^ < 3.0, and 0.01 < P^, = < 0.2, 

(A.29) 



where the critical properties of hydrocarbons are esti- 
mated by the following equations: 
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= 17.1419[exp( - 9.3145 x IQ-^TMeABP 
- 0.5444Sg + 6.4791 x IO'^'TmcabpSJ] 
X Ttif^irS^^'^K (A.30) 

= 4.6352 X 10*[exp( - 8.505 x IQ-^TMeABP 
- 4.8014S, + 5.749 x IQ-^rMeABpS^)] 



.4844_^4..0846 



(A.31) 



On the other hand, the viscosities of the components in 
the coke burning reaction mixture are expressed as func- 
tions of temperature only as follows 

= 6.476 x 10-« -I- 4.818 x IQ-^T, (A.32) 

^co = 5.38 x 10-* + 4.04 X 10-«T, (A.33) 

Hco, = 1.34 X 10"* + 4.5 X IQ-^T, (A.34) 

/iH^o = - 1.90 X IQ-^ + 3.85 X 10-«r, (A.35) 

yiN, = 1.101 X lO"' + 3.074 X 10"^r, (A.36) 

Finally, the mean viscosity of gaseous mixture is cal- 
culated by 

(A.37) 

The correlation equation for heat of vaporization of 
gas oil was obtained as Eq. (A.38) by nonlinear regression 
of literature data (Nelson, 1958). This equation is applic- 
able when the molecular weight of gas oil ranges from 
200 to 400. 



= 0.3843 



+ 1.0878X lO^exp - 



100 



5.153. (A.38) 



The vapor pressure of gas oil is calculated by solving 
for a trial-and-error solution of the following set of cor- 
relation equations adopted from literature (Technical 
Data Committee, 1988): 

P = 0 133322 X lQi^000.538!O-6.161S60)/{43ai-0.987672) 

for CO > 0.0022, 

Pg = 0.133322 X 10<26" l""-S.994296)/(95.-76«,-0.972546) 

for 0.0013 <cu < 0.0022, 

Pg = 0.133322 X 10<"-'«-085<o-6.412631)/(36<a-0.989679) 

for CO < 0.0013, (A.39) 



748.1 - 0.3 861 ri^.ABP 
^I&eABP = ^M«ABP "~ 1-3 889a 

x{Kf - 12)log(0.0098684PJ, 
a = 1 for TweABP > 477.8 K, 
a = 0 for Tmcabp < 333.3 K, 
a = (LSTmcabp - 659.7)/200, for 366.7 K 
^ Tmcabp < 477.8 K. 



(A.40) 



(A.41) 



(A.42) 



To calculate the three coefficients in the Antoine equa- 
tion for gas oil, vapor pressures are evaluated at three 
temperatures (TweABP, ^mcabp - 15, and Tmcabp + 15) 
using Eqs. (A.39)-(A.42) and are substituted into the 
following equations: 



FMi„B,g,Cig) = (Tw + C^)[^„ - log(P,.)] 
- B,g =0, j = 1,2,3. 



(A.43) 



The resulting set of equations are solved using the New- 
ton's method to get Aig, Big, and C^g. 

The thermal conductivity of hydrocarbons is needed to 
calculate thermal diffusivity and heat transfer coefficient. 
The following equation was obtained by nonlinear re- 
gression of the data in API technical data book (Technical 
Data Committee, 1988): 

fc^ = 1 X 10-^(1.9469 - 0.374Mh,„ 

+ 1.4815 X 10" + 0.1028r). (A.44) 

Then the interface heat transfer coefficient between solids 
and hydrocarbon gases can be estimated from the cor- 
relation between the Nusselt and Reynolds numbers 
(Kunii &. Levenspiel, 1991): 



Then the overall interchange coefficient Kj for mass 
transfer between the emulsion and the bubble phases can 
be computed using the interchange coefficient Ki be- 
tween bubble and cloud and between cloud and 
emulsion (Kunii & Levenspiel, 1991) 
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K,=4.5^ + 5.85[-^j (A.47) 

where the effective diffusion coefficient of gaseous 
mixture in the reactor and regenerator is calculated using 
the correlation equation proposed by Baird and Rice 
(1975): 

Df = 0.35(guJ''3Z)'*'^ {A.49) 



Appendix B. Additional simulation data 

The following is a list of additional data used in 
our simulation. These data together with those given 
in Tables 1-4 correspond to the base case simulation 
condition. 

/. Initial conditions: 

C,^j,r = 0.009, w,«r = 485.0 kg, 

w,j,r = 38 000.0 kg, y„„« j- = 0.27, 

y,,RT = 0.51, y,,«r = 0.14, T^r = 787.0 K, 

w,^G - 990 kg, w,„cj = 182 000 kg, 

C,,B = 0.001, Tj,= 991.0 K, 

Co,£ = 0.0005 kg mol/m^ Ccoe = 0.0003 kg mol/m^ 
Cco^E = 0.004 kg mol/m^ Ch.oe = 0.003 kg mol/m^ 
Cnz£ = 0.02 kg mol/m^ Cq.b = 0.0005 kg mol/m^ 
CcoB = 0.0003 kg mol/m^, Cco.b = 0.004 kg mol/m^ 
Ch.ob = 0.003 kg mol/m^, Cn^b = 0.02 kg mol/m^ 
Co,F = 0.0001 kg mol/m^, Ccof = 0.0002 kg mol/m^, 
Cco.F = 0.004 kg mol/m^, Ch.of = 0.003 kg mol/m^ 
CN,f = 0.02 kg mol/m^ Tp = 1010.0 K, 
C,;t^. = 0.0008. 



2. Fia/ue parameters: 
k„MF = 4.35 kg/(skPa<^-'), 
fc„cY2 = 140.0 kg/{sm^-'), 
=69.0kg/(skPa<'-'), 
a„Mr = 1-0, a„cYi = 1-0, 
a.,svi = 0.9, a„sK2 = 0.9. 



fe„cYi = 120.0 kg/(sm^-^), 
= 4.25 k^(s kPa°-^), 
Ksvi =69.0kg/(skPa°-^), 
a„cY2 = 1.0, a,sG = 0.8, 



5. Heat transfer coefficients: 

U^r = 0.035 kJAm^ s K), Uo = 0.05 kJAm^ s K), 
Up = 0.05 kJ/(m' s K). 

4. Miscellaneous parameters: 

CckSTO = 0.0004, fe,,o = 0.0054, = 0.3311 s/kg, 

AH,,,i^sT = 0 kJ/kg, PfcCYi = 350 kg/s, 
P6CY2 = 350 kg/s, Picii = 560 kg/s, 
^6cz,2 = 560 kg/s, cr = OK, Bg„ = 0.6, 

= 1.5 s-^ vcYi = 1. 
rjcyi = 1, AT.cYi = 10, Ar,cY2 = 10. 
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(7-68) 



u.= Ufd- 0.06384 Vp') 
Co = 0.44 

7-12] W»m?1-l=*^J .-^+»^^^^^° 



(7-69) 





33370 
11870 
6220 
4900 
412 
344300 
401072 
56772 



2090 
8960 
19040 
10350 

7950 = 6>>«-'- '7?o 

4900 
412 

347370 5 34V^.-.+ lo 
401072 
53702 



94.265 
23.98 
18.40 
272.22 
14.207 

423.072 
423.072 



104.657 
182.22 
184.854 
50.735 
23.521 
272.22 
14.207 

832.414 
832.414 



^^|^'AP = 41 .4 + 2.5x30x 10"' =42.15kPa 

JPi55^f!;m54-|S[^J»JH£j^ 4 X 560/100 = 22. 40kPa 

1000 + 19 . 75 - 253 . 75kPa = 0 . 254MPa 



0.234 X 



22.4 X [(579 + 273)/273] x [ 1 .033/(0.254 x 10.19)] 
= 11804.07m'/h = 3.28mVs 



{i^ =423. 012 > 



^^mAam^^^ u^= =3.28/0.567 = 5. 78m/so 

=832.414 X 22.4 X [(537 + 273)7273] x [l .033/0.24 x 10.19] 
= 23368 . 17mVh = 6 . 49mVs 

V^/A =6.49/0.567 = 11. 45m/s 
5)M^^¥i^^i^o 



^"T78 

6)-?S^^^o 

r tt] ^ S » m S ^' yS ^ eTl * 6^ b t ttlli ;^ V* ^ ^ S o 

ia P , = 56772/1 1804 . 07 = 4 . 8 lkg/m\ 
J^^U,p^= 53702/23368 .117 = 2. 30kg/tn' , 

2.M?t-Wm»i+» 
(l)tt^@MMS 









^±5 




5.78 


11.45 


8.29 




2.5 


1.5 


1.86 




2.31 


7.63 


4.45 



= 1.3 X 10-' (30/0. 85) X 3.4x8.29" 
= ,107.21Pa 

A/*^ = X I X g = 3 .4 X 30 X 9. 81 = 1000 .6Pa 
AP^ =0.5 xpg xV = 0.5 X 3.4x8. 29" = 14d_4Pa- 
AP^ = G3 

G3 = 344300/(0.567 x 3600) = 168.7kg/(m' -s) 

AP^ = 168.7 X 7.63 = 1286. 99Pa 

AP^ = 5.7x 10-' X C.Lxgx (gdj)-°' 

= 5.7x 10"' X 168.7x30x9.81 X (9.81 x 0.85)"° ' 

= 980.02Pa 

AF^ + AP^ + AP^ + AP^ + AP^f = 3521 .24Pa 

242 



ANNEX IV 

Process Calculation and Technical Analysis of 
Catalytic Cracking 

Editor-in-chief: Cao Hanchang Hao Xiren Zhang Han 



Petroleum Industry Press 



1 



(6) Gas velocity head at the outlet is represented by AP^^ and calculated 
according to the following equation: 

2. Pressure drop of horizontally transported dilute phase 

There are a lot of equations for calculating pressure drop of horizontally 
transported dilute phase, of which the one available hereinbelow is 
recommended: 

1 2 ^ 2/,p,u?L r, .y>L^1 (7-69) 

J' S Pp d,p u. 
u. = ur(l-0.0638<Vr) 
Cd = 0.44 

In addition, when f^ujif^u^ > 1, the value of 1 is taken for use in the 
above equation. 

IV. Accounting of the Operational Parameters of Riser reactors 

[Example 7-12] The operational parameters are calculated based on a 
sample riser. 

Original Data 



1 . Reaction Conditions 



Temperature at the outlet of the riser 


537 °C 


Temperature of the regenerated catalyst 


686 °C 


Preheating temperature of the fresh feedstock 


256 °C 


Temperature of the recycle slurry 


360 °C 


Temperature of the heavy cycle oil(HCO) 


325 °C 


Pressure at the top of the disengager 


0.234 MPa (absolute 




pressure) 
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Xot3.1 flow rate of stcarti in the riser 


4900 kg/h 


Flow ra.te of the distillate oil 


33370 kg/h, 94.265 kmol/h 


i/iuw I aLC Lii Liic ucav y on 


11870 kg/h, 23.98 kmol/h 


-TiLiw xaLC oi Lixc icoyv./ic &iuiiy 


6220 kg/h, 18.4 kmol/h 


Flow rate of the heavy cycle oil(HCO) 


0 


Circulation rate of the catalyst 


344300 kg/h 


Flow rate of the flue gas carried by the regenerated 


412 kg/h, 14.207 kmoVh 


catalyst 




2. Flow Rates of the Products 


Dry gas and loss 


2090 kg/h, 104.657 kmol/h 


LPG 


8960 kg/h, 182.22 kmol/h 


Gasoline 


19040 kg/h, 184.854 kmol/h 


Light cycle oil 


10350 kg/h, 50.735 kmol/h 


Slurry 


1730 kg/h 


Coke 


3070 kg/h 


3. Properties of the Feedstock and Products (omitted) 


4. Dimensions of the Riser Structure 


Inner diameter 


0.85 m 


Length of the prelift zone 


4.0 m 


Length of the horizontal zone 


5.0 m 


Length of the vertical zone 


30 m 


5. Parameters of the Riser such as Pressure Drop 


Total pressure drop (as actually measured) 


41.4 kPa 


Pressure drop of Horizontal Section and cyclone as 


5.7 kPa 


actually measured 




Density of the prelift zone 


560 kg/m' 


Density of the counter draught 


2.5 kg/m' 


Dilute-phase static pressure from the top of the 


0.5 kPa 


disengager to the outlet of the riser 





1 . Material Balance in the Riser and Calculation of Velocity and Density 
1) Material Balance 
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Materials at the Inlet of the Riser 



Name 


Ajfrt**,*^ T71^,-., /l>-#*/U\ 

Mass rlOW {Kg/n) 


Molar Flow (kmol/h) 


Distillate oil 


33370 


94 265 


Heavy oil 


11870 


23.98 


Heavy cycle oil/ recycle slurry 


6220 


18.40 


steam 


4900 


272.22 


Flue gas carried by the 


412 


14.207 


regenerated catalyst 






Catalyst 


344300 




Total 


401072 


423.072 


Total of vapor 


56772 


423.072 


Materials at the outlet of the riser 


Name 


Mass Flow (kg/h) 


Molar Flow (kmol/h) 


Dry gas and loss 


2090 


104.657 


LPG 


8960 


182.22 


Gasoline 


19040 


184.854 


Light cycle oil 


10350 


50.735 


Slurry 


7950 


23.521 


Steam 


4900 


272.22 


Flue gas carried by the 


412 


14.207 


regenerated catalyst 






Coke/Catalyst 


347370 




Total 


401072 


832.414 


Total of vapor 


53702 


832.414 



2) Pressure and Temperature at the Inlet of the Riser 
(1) Pressure 

Operational pressure at the inlet of the riser = Pressure in the disengager 
+ Dilute-phase static pressure from the top of the disengager to the outlet 
of the riser + Total pressure drop in the riser (excluding the prelift zone) 
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The measured value of the total pressure drop in the riser is 41.4 kPa, 
which is revised as: 

AP = 41 .4 + 2.5 X 30 X 10"^ = 42.15 kPa. 

The static pressure in the prelift zone that should be deducted is 4 x 
560/100 = 22.40 kPa. 

Then, the total pressure drop in the riser excluding the prelift zone is 
19.75 kPa. Consequently, the operational pressure at the inlet of the riser 
is: 

0.234 X 1000 + 19.75 = 253.75 kPa = 0.254 MPa. 
(2) Temperature 

Owing to heat balance, the mixed temperature obtained after gasification 
of the feedstock resulting from the catalyst-oil contact is 597 °C (the 
calculation procedure is omitted here). 

3) Linear Velocity at the Inlet of the Riser 

The volume flow rate of the vapor at the inlet of the riser is calculated as 
follows: 

V^, = 423.072 X 22.4 x [(597 + 273) / 273] x [1.033 / (0.254 x 10.19)] = 
11804.07 mVh = 3.28 mVs. 

The linear velocity at the inlet of the riser is: u^- = V^J A = 3.28/0.567 = 
5.78 m/s. 

4) Linear Velocity at the Outlet of the Riser 

The volume flow rate of the vapor at the outlet of the riser is calculated as 
follows: 

V^, = 832.414 X 22.4 x [(537 + 273) / 273] x [1.033 / 0.24 x 10.19] = 
23368.17 mVh = 6.49 mVs. 



Thus, the linear velocity at the outlet of the riser, u^^, is calculated as: 
w.e =V,JA = 6.49/0.567 = 1 1 .45 m/s. 



5) Mean Linear Velocity in the Riser 



It is calculated based on the logarithmic mean value: 

i 11.45-5.78 



-. = --ir^= , 11.45 =8-29-^^ 



6) Density of the Vapor 



The ratio of the weight flow rate of the vapor to the volume flow rate of 
the vapor is the density of the vapor. 

With respect to the inlet: Pgj = 56772/11804.07 = 4.81 kg/m^; 
With respect to the outlet: Pg,= 53702/23368.117 = 2.30 kg/m\ 

The logarithmic mean density in the riser is: Pg = 3.4 kg/m^. 

2. Calculation of the Inventory of the Riser 

(1) Estimation of the Velocity of Solid Particles 



The velocity of the particles is obtained by estimating the sliding 
coefficients at the inlet and outlet in light of relevant information and the 
gas velocity in the riser. 





Inlet 


Outlet 


Average 


Velocity of the vapor, m/s 


5.78 


11.45 


8.29 


Estimated K^i 


2.5 


1.5 


1.86 


Velocity of the catalyst, m/s 


2.31 


7.63 


4.45 



In this step, only the catalyst-accelerated pressure drop is calculated, and 
the results need to be revised further. 



(2) Calculation of Pressxire Drops Excluding the Catalyst Static Head 
APgf = 1 .3 X 10-2 ^^/^^^ ^ ^2 

= 1.3 X 10-2 (30/0.85) X 3.4 X8.292 

= 107.21 Pa 

APgh = PgXLxg = 3.4x30x9.81 - 1000.6 Pa 
APg, = 0.5 X pg X ul =0.5 X 3.4 X 8.29^ = 146.4 Pa 

G, = 344300/(0.567 x 3600) = 168.7 kg/(m2.s) 

AP,, = 168.7 X 7.63 = 1286. 99 Pa 

AP,f = 5.7 X 10-2 X G/. X g X (g^^)-o.5 

= 5.7 X 10-2 X 168.7 X 30 X 9.81 x (9.81 x 0.85)-^ ^ 
= 980.02 Pa 

Thus, APgf + APgh + APg, + AP,, + AP,f = 3521 .24 Pa 



